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CHAPTER  I 


INTRODUCTION 


Recently  at  the  Applied  Research  Laboratory  of  The  Pennsylvania 
State  University  (ARL/PSU) ,  a  cantilever-type  balance  was  designed  to 
measure  the  lift,  drag,  and  pitching  moment  characteristics  of  a  hydro¬ 
foil  in  cavitating  flow.  As  a  preliminary  verification  of  the  test 
setup,  several  standard  NACA  (National  Advisory  Committee  for 
Aeronautics)  airfoils  were  tested.  Data  from  these  tests  were  compared 
to  those  previously  published  by  NACA  and  by  those  measured  by 
Kermeen  (l)  at  the  California  Institute  of  Technology.  The  ARL/PSU  test 
results  differed  substantially  from  both  sets  of  previously  published 
data.  Kermeen  did  not  specifically  address  the  problem  of  water  tunnel 
balance  design  because  he  used  a  balance  of  established  accuracy. 
However,  there  are  also  slight  differences  between  the  results  Kermeen 
reported  and  those  reported  by  NACA.  Because  of  these  considerations,  a 
program  was  initiated  to  study  the  problem  of  hydrofoil  balance  design 
with  a  particular  emphasis  on  the  action  of  the  tunnel  walls  and  viscous 
effects  on  the  measured  forces. 

The  ultimate  goal  of  the  program  was  to  develop  the  capability  of 
accurately  measuring  force  coefficients  characteristic  of  new  hydrofoil 
profiles  in  fully  wetted  and  cavitating  flows.  To  achieve  this  objec¬ 
tive,  five  intermediate  objectives  were  established. 

The  first  of  these  objectives  was  to  determine  why  the  12-in. 
water  tunnel  balance  failed  to  measure  the  forces  accurately. 
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The  second  objective  was  to  build  and  test  a  balance  which  was 
similar  in  principle  to  the  ARL/PSU  water  tunnel  balance  for  use  in  the 
48-in. -diameter  wind  tunnel.  By  use  of  the  48-in.  wind  tunnel,  larger 
airfoil  geometries  could  be  employed  in  an  effort  to  duplicate  water 
tunnel  Reynolds  numbers  (Re).  Waterproofing  problems  and  the  large 
forces  involved  in  water  flows  are  also  eliminated  in  wind  tunnel  test¬ 
ing. 

The  third  objective  was  to  determine  what  factors  (and  their  magni¬ 
tude)  affect  flow  past  a  hydrofoil  in  a  two-dimensional  channel.  The 
understanding  of  factors  affecting  the  flow  can  possibly  give  an  insight 
as  to  why,  hydrodynamically  speaking,  the  original  ARL/PSU  balance 
failed.  Such  knowledge  can  also  provide  a  sound  information  base  for 
designing  or  modifying  hydrofoil  balances. 

The  fourth  objective  was  to  determine  a  theoretical  and/or  empiri¬ 
cal  correction  procedure  to  use  with  the  existing  balance,  one 
accounting  for  all  important  physical  interactions  in  the  channel  and 
converting  the  balance  measurements  into  accurate  two-dimensional 
forces.  Depending  on  the  adequacy  of  the  correction  procedure,  a 
decision  would  be  made  whether  to  continue  to  use  the  existing  balance 
in  conjunction  with  a  correction  procedure  or  to  redesign  the  balance 
based  on  new  knowledge. 

The  fifth  objective  was  to  formulate  a  standard  procedure  for 
measuring  and  reducing  data  with  the  ultimate  tunnel-balance 
configuration  at  ARL/PSU.  This  procedure  would  be  documented  to  provide 
a  reference  for  future  two-dimensional  hydrofoil  testing. 


CHAPTER  II 


BACKGROUND 

Although  air  and  water  are  both  fluids,  there  are  fundamental  dif¬ 
ferences  between  air  and  water  tunnel  testing.  Consideration  of  these 
differences  led  to  the  conclusion  that  a  cantilever  balance  system  as 
described  in  Chapter  III  is  the  best  choice  for  water  tunnel  testing. 

A  traditional  method  of  measuring  force  coefficients  in  wind  tun¬ 
nels  has  been  by  use  of  pressure  taps.  In  a  two-dimensional  channel, 
lift  being  created  by  an  airfoil  exerts  a  force  on  the  ceiling  and  floor 
of  the  tunnel.  By  placing  pressure  taps  on  these  walls  and  integrating 
the  pressures,  the  lift  can  be  determined.  Drag  forces  can  also  be 
determined  by  placing  pressure  taps  along  the  surface  of  the  airfoil  and 
integrating  the  readings  in  the  drag  direction.  Another  method  of  meas¬ 
uring  drag  with  pressure  sensitive  instruments  is  to  measure  the 
momentum  deficit  created  by  the  airfoil  in  a  two-dimensional  channel. 
This  measurement  is  accomplished  by  traversing  a  pressure  probe  through 
the  airfoil  wake  and  integrating  the  readings. 

Testing  airfoil  shapes  in  water  (hydrofoils)  introduces  additional 
problems.  Among  these  problems  are  handling  larger  gross  forces  and 
waterproofing  requirements.  The  major  area  of  concern,  however, 
involves  cavitation.  When  hydrofoils  are  tested  in  water  tunnels,  they 
may  operate  in  three  different  flow  regimes.  The  first  regime  is 
noncavitating  or  fully  wetted  flow.  This  range  of  flow  is  that  which 
falls  within  the  scope  of  the  present  study,  and  it  differs  little  from 
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low-speed  wind  tunnel  testing.  The  second  regime  is  the  partial  cavita¬ 
tion  range  or  that  regime  where  the  cavitation  bubble  begins  and  ends  on 
the  surface  of  the  body.  The  third  regime  is  fully  cavitating  flow  or 
cavity  flow.  The  present  research  concentrates  on  the  first  flow 
regime . 

A  balance  system  and  correction  procedure  for  water  tunnel  testing 
must  be  capable  of  accurately  measuring  forces  in  all  three  flow 
regimes.  Compared  to  wall  effects  in  fully  wetted  flows,  the  wall 
effects. on  fully  cavitating  profile  forces  are  not  as  pronounced  pro¬ 
vided  that  the  cavitation  number  is  based  on  measured  cavity  pressure. 
Elaborate  cavity  flow  calculations  can  be  used  to  correct  the  force  data 
in  this  last  case.  Although  sidewall  boundary  layer  effects  are  still 
observed  in  the  fully  cavitating  flow  regime,  past  experience  indicates 
that  they  have  a  negligible  effect  on  the  measured  forces  in  this  case. 
The  effects  of  tunnel  walls  in  the  partially  cavitating  regime  are 
beyond  the  scope  of  this  study. 

Placing  pressure  taps  on  a  hydrofoil’s  surface  or  on  the  tunnel 
walls  will  cause  premature  cavitation  and  cavitating  taps  have  been 
shown  to  give  spurious  pressure  readings.  This  is  obviously  unaccept¬ 
able,  especially  when  tests  are  being  conducted  in  the  partial 
cavitation  regime.  For  these  reasons,  water  tunnel  testing  is  best 
performed  by  directly  measuring  the  forces  with  a  balance. 

With  a  mechanical  balance,  aerodynamic  forces  can  be  measured 
directly  without  pressure  taps  marring  the  model  surface.  Mechanical 
balances  are,  however,  not  totally  free  of  problems.  Balances  measure 
all  forces  applied  to  a  model.  If  forces  occur  on  a  model  which  are  not 
those  associated  with  two-dimensional  flow  (due  to  such  causes  as 


secondary  flows  where  the  model  intersects  a  tunnel  wall),  the  balance 
will  also  measure  them.  A  correction  factor  is  then  required  to  reduce 
this  total  force  to  a  two-dimensional  force.  Despite  the  required  cor¬ 
rections,  mechanical  balances  appear  to  be  the  best  method  for  measuring 
forces  on  cavitating  hydrofoils.  Previous  work  by  various  authors  sub¬ 
stantiates  this  point  of  view. 

2.1  Baseline  Data 

As  a  prelude  to  the  present  study,  it  would  be  beneficial  to  exam¬ 
ine  published  wind  tunnel  data  to  establish  a  credible  reference  frame 
for  the  data  of  this  study.  A  wealth  of  airfoil  test  data  were  measured 
by  NACA  at  the  Langley  and  Ames  research  centers  in  the  1930's  and 
1940's.  Loftin  and  Smith  (2)  conducted  numerous  tests  on  a  variety  of 
two-dimensional  airfoil  shapes  at  Langley  Aeronautical  Laboratory  in 
1949.  They  used  airfoils  of  3-ft  chord  which  completely  spanned  a 
rectangular  test  section  of  3-ft  width  and  7. 3-ft  height.  Lift 
measurements  were  made  by  taking  the  difference  between  integrated 
pressure  reactions  on  the  floor  and  ceiling  of  the  tunnel.  Lift  was 
also  measured  by  use  of  a  three-component  force  balance.  Small  gaps 
were  allowed  between  the  walls  and  airfoil  tips  to  permit  freedom  of 
movement  of  the  model.  Although  these  end  gaps  are  a  source  of  error, 
the  two  sets  of  lift  measurements  agreed  within  the  experimental  error 
of  the  test.  Drag  measurements  were  made  with  a  wake  survey  apparatus. 
The  model  end  gaps  were  sealed  with  felt  packing  which  reduced  errors 
due  to  induced  drag.  Pitching  moment  was  measured  by  a  torsional 
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Jacobs  and  Sherman  (3)  tested  numerous  airfoil  profiles  as  a  func¬ 
tion  of  Reynolds  number  (Re) .  They  tested  at  Re  values  as  iow  as 
170,000.  These  tests  were  within  the  range  of  Re  achievable  in  the 
48-in.  wind  tunnel  at  ARL/PSU.  The  airfoils  tested  in  Reference  (3) 
were  finite  aspect  ratio  models  which  did  not  completely  span  the  wind 
tunnel.  In  effect,  three-dimensional  data  were  measured  and  corrected 
to  two-dimensional  data.  The  authors  reported  that  the  magnitude  of 
experimental  errors  increased  as  Re  was  lowered.  They  further  stated 
that  drag  and  pitching  moment  results  became  relatively  inaccurate  below 
Re  =  800,000  due  to  limitations  imposed  by  the  sensitivity  of  the  meas¬ 
uring  equipment.  Inaccuracy  of  these  quantities  became  so  pronounced 
that  airfoil  characteristics  dependent  on  drag  and  pitching  moment 
results  (optimum  lift  coefficient,  aerodynamic  center  of  pressure,  etc.) 
were  considered  unreliable  and  in  most  cases  were  not  presented  below 
Re  =  800,000.  The  data  presented  in  Reference  (3)  should  be  verified  as 
an  integral  step  in  this  study. 

2.2  Previous  Water  Tunnel  Results 

Kermeen  (l)  conducted  a  series  of  experiments  to  determine  whether 
accurate  force  measurements  could  be  obtained  with  a  recently  modified 
hydrofoil  force  balance  (4)  at  the  California  Institute  of  Technology 
(CIT) .  He  used  a  cantilever  hydrofoil  balance  arrangement  (Figure  1) 
and  tested  NACA  4412  and  Walchner  profile  7  hydrofoils. 

The  airfoil  was  originally  mounted  on  a  stub  spindle  as  depicted  on 
the  left  of  Figure  2.  This  arrangement  necessitated  small  clearance 
gaps  on  each  end  of  the  airfoil.  Kermeen  stated  that  this  arrangement 
caused  problems.  There  were  radial  interferences  due  to  deflection  of 


NACA  4412  hydrofoil  models  used  by  Kermeen.  The 
original  stub  spindle  is  shown  on  the  left.  The 
newer  spindle-disk  assemply  is  shown  on  the  right. 


the  two-dimensional  channel  walls.  There  was  also  a  problem  of  water 
jetting  through  the  stub  spindle  clearance  gap  from  the  high-pressure 
dead-water  region.  For  these  reasons,  Kermeen  switched  to  the  enlarged 
spindle  disk  arrangement  pictured  on  the  right  of  Figure  2.  The  5-in.- 
diameter  disk  was  attached  to  a  hydrofoil  of  3-in.  chord,  thereby  elim¬ 
inating  the  end  gap  on  one  side  of  the  hydrofoil. 

This  new  arrangement  eliminated  corrections  that  may  have  been 
necessary  due  to  the  end  gap  on  the  balance  side  of  the  model.  Kermeen 
gives  a  good  description  of  five  of  the  traditional  sources  of  error  in 
wind  and  water  tunnel  testing: 

1.  Solid  blockage.  This  considers  the  effect  of  the  model 
physically  blocking  the  flow  and  subsequently  changing  the 
axial  velocity  past  the  hydrofoil. 

2.  Wake  blockage.  This  is  a  similar  effect  and  accounts  for 
the  reduction  of  velocity  inside  the  hydrofoil  wake. 

3.  Lift  effect.  This  accounts  for  the  compression  of  the 
streamlines  around  the  hydrofoil  due  to  the  tunnel  walls. 

The  effect  is  treated  as  a  change  in  hydrofoil  camber  and 
incidence  angle. 

4.  Horizontal  buoyancy.  This  correction  compensates  for  the 
additional  drag  created  on  the  model  due  to  the  existing 
static-pressure  gradient  in  a  two-dimensional  channel. 

5.  Boundary  layer  interference.  This  accounts  for  the  depar¬ 
ture  from  two-dimensional  flow  near  the  walls  due  to  the 
interaction  of  the  wall  boundary  layer. 

Kermeen  also  considered  the  effect  of  skin  friction  on  the  disk. 

He  mounted  the  foil  on  the  opposite  side  of  the  channel  and,  with  a 


small  clearance  between  the  disk  and  foil  (approximately  0.002  in.), 
measured  the  forces  sensed  by  the  disk  alone  for  various  angles  of 
attack.  These  tare  forces  were  then  subtracted  from  the  appropriate 
test  data. 

In  addition  to  the  above  consideration,  Kermeen  further  examined 
the  effect  of  flow  through  the  gap  between  the  hydrofoil  end  and  the 
tunnel  wall.  With  a  gap  as  small  as  0.002  in.,  the  author  reported  that 
flow  through  this  gap  was  clearly  visible.  He  varied  the  gaps  from 
0.001  in.  to  0.032  in.  and  for  an  angle  of  attack  of  4°  found  that  the 
drag  coefficient  increased  by  25%.  The  lift  was  reduced  by  8%  for  these 
same  conditions.  The  moment  coefficient  appeared  to  be  unaffected  by 
gap  size. 

In  his  final  data  reduction,  Kermeen  chose  to  include  corrections 
for  all  of  the  tunnel  interference  effects  except  boundary  layer  inter¬ 
ference.  He  felt  that  the  boundary  layer's  effect  on  lift  was  negli¬ 
gible  and  he  made  no  statement  about  its  effect  on  drag.  He  also 
included  corrections  for  the  tare  force?  on  the  disk.  He  did  not  spe¬ 
cifically  state  whether  or  not  he  included  corrections  for  the  gap 
effect.  Kermeen's  lift  and  moment  data  for  the  NACA  4412  profile  agreed 
very  closely  to  the  NACA  data  (3) .  The  drag  data  agreed  well  in  the 
lower  drag  range  but  for  attack  angles  greater  than  5°  Kermeen's 
results  showed  higher  drag  values.  The  thrust  of  Kermeen's  study  was  to 
demonstrate  as  clearly  as  possible  that  the  balance  could  accurately 
measure  hydrodynamic  forces.  This  enabled  him  to  compare  his  new  data 
on  cavitating  and  noncavitating  Walchner  hydrofoils  with  Walchner's 
results,  which  were  not  very  accurate.  He  did  not  attempt  to  determine 
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what  factors  beyond  the  classical  corrections  affected  balance  perform¬ 
ance. 

Earlier  studies  of  a  similar  nature  were  conducted  by  Daily  (5)  at 
the  California  Institute  of  Technology  (CIT) .  He  ran  a  series  of  tests 
on  a  NACA  4412  hydrofoil  using  a  three-component  force  balance  (from 
which  the  balance  of  Kermeen's  report  was  derived)  in  the  high-speed 
water  tunnel. 

Daily's  setup  is  depicted  in  Figure  3.  He  used  a  cantilever 
arrangement  with  a  small  supporting  spindle.  The  size  of  the  supporting 
spindle  resulted  in  clearance  gaps  on  each  end  of  the  model.  Because  of 
the  magnitude  of  the  forces  at  the  higher  angles  of  attack,  Daily  used  a 
full-span  and  a  split-span  hydrofoil.  Splitting  the  span  in  half,  as 
indicated  in  Figure  3,  enabled  the  balance  to  sense  only  half  of  the 
gross  forces  on  the  hydrofoil.  The  full-span  model  was  used  for  testing 
at  the  lower  angles  of  attack  and  the  semi-span  version  was  used  for 
testing  in  the  higher  lift  range. 

Daily  implied  that  his  data  required  few  corrections.  He  reported 
that  boundary  layer  effects  in  the  region  where  the  hydrofoil  meets  the 
tunnel  wall  should  be  neglected.  Since  the  foil  was  mounted  only  one 
tunnel  diameter  from  the  final  contraction  of  the  flow,  the  boundary 
layer  is  relatively  thin  when  it  interacts  with  the  foil.  He  further 
reported  that  incipient  cavitation  photographs  suggested  a  uniform  span- 
wise  flow.  Daily  also  claimed  that  the  clearance  gap  between  the  tunnel 
walls  and  the  model  test  span  was  small  enough  to  make  its  effect  negli¬ 
gible.  He  further  stated  that  wall  interference  or  "blocking"  was 
negligible  at  low  angles  of  attack  and  increased  with  increasing  angles. 
The  magnitude  of  the  blocking  effect  was  not  evaluated  for  his  tests. 
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FULL  SPAN  HYDROFOIL  SPLIT  SPAN  HYDROFOIL 


Figure  3.  Sketches  Illustrating  Daily's  balance-hydrofoil  system. 
The  top  sketch  shows  a  schematic  diagram  of  the 
balance.  The  bottom  sketch  shows  details  of  the  model 
installations  in  the  tunnel  working  section. 
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Daily  compared  his  test  results  to  the  NACA  wind  tunnel  tests  of 
Reference  (3) .  He  generally  found  good  agreement  at  low  angles  of 
attack  but  his  measurements  differed  by  as  much  as  10%  at  the  higher 
angles  of  attack.  He  further  concluded  that  the  low  lift  range  was  the 
operating  range  of  many  hydrofoil  applications  and  that  with  his  test 
setup  good  comparative  results  could  be  obtained. 

2.3  Recent  ARL  Water  Tunnel  Tests 

In  the  previously  mentioned  test  program  at  ARL/PSU,  the  work  of 
Kermeen  and  Daily  was  used  as  a  guide  for  designing  the  hydrofoil  and 
balance  hardware.  The  ARL/PSU  balance  and  a  typical  hydrofoil  are  shown 
in  Figure  4.  This  three-component  balance  measured  lift,  drag,  and 
pitching  moment  by  use  of  strain  gages  attached  to  tension  members  which 
are  aligned  axially  and  normally  to  the  hydrofoil  (Figure  4) .  The  spin¬ 
dle  attachment  on  the  foil  was  a  compromise  between  Kermeen's  larger- 
than-chord  disk  and  Daily’s  equal-to-hydrofoil-thickness  spindle.  The 
aspect  ratio  (AR)  of  the  ARL/PSU  hydrofoils  was  0.9  as  compared  to 
Daily's  AR  *  3.33  and  Kermeen's  AR  “  0.97. 

In  the  original  ARL/PSU  program,  a  NACA  4412  hydrofoil  was  tested 

(Figures  5A  and  3B) .  The  pitching  moment  results  agreed  extremely  well 

with  published  NACA  wind  tunnel  data  (3) .  The  slope  of  the  lift  curve 

and  CT  were  approximately  20%  lower  than  published  data.  The  drag 
max 

curve  was  within  25%  of  published  data  at  the  zero-lift  condition. 
However,  at  all  other  values  of  lift,  the  corresponding  drag  values  were 
higher  than  published  data  by  as  much  as  a  factor  of  ten.  These  conclu¬ 
sions  were  obtained  from  uncorrected  data  taken  just  prior  to  the 
termination  of  the  force  phase  of  the  test  program. 
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Several  attempts  were  made  to  correct  the  problems  with  the  bal¬ 
ance.  When  the  balance  was  mounted  in  the  horizontal  configuration,  it 
was  noted  that  the  hydrofoil  deflected  severely  under  loads  due  to  weak¬ 
nesses  in  the  flexures  between  the  force  cubes.  In  order  to  eliminate 
this  problem,  the  balance  was  rotated  90°  so  that  the  cubes  were  in  a 
vertical  configuration.  This  rotation  eliminated  the  model  deflection 
problem  but  required  the  drag  or  axial  force  to  be  measured  as  the  sum 
of  two  cube  readings  instead  of  one.  At  high  angles  of  attack,  the 
axial  force  was  small  compared  to  the  normal  force.  This  vertical 
arrangement  meant  that  a  small  force  was  being  divided  in  half  and  that 
drag  was  being  measured  as  the  sum  (or  difference,  depending  on  the 
signs  of  the  calibration  constants)  of  two  small  values.  The  realign¬ 
ment  of  the  balance  was  one  of  several  attempts,  all  of  which  failed  to 
solve  the  problem. 

2.4  Further  Pertinent  Data 

More  recently,  an  extensive  investigation  of  profile  forces  was 

conducted  by  Otsuka  and  Sugiyama  [6) .  They  conducted  numerous  tests  on 

an  airfoil  spanning  a  rectangular  channel  with  a  variable  clearance 

between  the  airfoil  tip  and  the  channel  wall.  They  measured  section 

lift  and  drag  forces  at  various  spanwise  locations  while  varying  the 

airfoil  angle  of  attack,  boundary  layer  thickness,  and  tip  clearance. 

They  placed  numerous  pressure  taps  on  the  airfoil  surface  and  integrated 

pressure  readings  to  determine  the  forces.  An  interesting  feature  of 

their  data  (on  a  plot  of  C^-versus-span  location)  was  that  accurate 

section  C,  values  could  apparently  be  obtained  for  a  given  airfoil  when 
a 

the  measuring  probe  was  far  enough  away  from  the  tunnel  wall.  This 
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is  the  fact  which  earlier  NACA  investigators  had  used  in  order  to  obtain 
reliable  airfoil  section  data. 

2.5  Discussion  of  Data  Base 

Consideration  of  the  work  of  these  previous  authors  yielded  several 
salient  points  which  were  useful  in  designing  a  program  to  accomplish 
the  goals  of  this  report. 

First,  in  order  to  establish  correction  factors  for  a  given  balance 
design,  there  must  exist  some  absolute  data  base  for  comparison  pur¬ 
poses.  For  a  given  airfoil,  the  lift,  drag,  and  pitching  moment  charac¬ 
teristics  must  be  known  before  the  magnitude  of  the  balance  corrections 
can  be  determined.  The  NACA  data  for  the  Reynolds  numbers  being  consid¬ 
ered  here  (generally  Re  <_  1,000,000)  appeared  to  be  insufficient. 
Therefore,  it  was  necessary  to  generate  a  valid  data  base  for  comparison 
purposes. 

The  second  point  concerned  the  approach  used  by  NACA  in  obtaining 
their  data.  Lift,  drag,  and  pitching  moment  were,  in  most  cases, 
measured  independently  and  by  the  most  accurate  technique  for  each 
force.  Lift  and  pitching  moment  could  be  measured  accurately  by  force 
balances,  while  drag  could  be  more  accurately  obtained  by  momentum 
surveys  via  pressure  surveys.  The  idea  of  measuring  quantities 
separately  (using  the  most  accurate  technique  for  each)  and  combining 
data  seems  to  be  an  excellent  solution  to  the  problem  of  establishing  a 
reference  data  base.  Data  shown  in  Figure  6,  taken  from  Kermeen, 
supported  this  idea.  Lift  and  pitching  moment  coefficients  for  a 
NACA  4412  airfoil  are  presented.  Consideration  must  be  given  to  the 
spread  of  Reynolds  numbers.  However,  all  factors  considered,  the  data 
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are  in  excellent  agreement.  Since  the  different  experimental  results 
agree  well  with  each  other  and  Reference  (2)  states  that  lift  measured 
by  both  pressure  surveys  and  a  force  balance  shows  no  significant  dif¬ 
ference,  the  logical  conclusion  is  that  the  lift  values  for  a  set  of 
reference  data  can  be  obtained  by  measuring  with  a  force  balance. 

Roughly  the  same  argument  also  holds  for  measuring  pitching  moment  coef¬ 
ficients. 

The  drag  measurements  required  special  consideration.  The  polars 
of  Figure  7  show  that  there  is  poor  agreement  in  measured  drag  coef¬ 
ficients  between  the  several  experiments.  Because  of  the  small  forces 
involved  and  the  sensitivity  of  drag  to  various  factors  (gap  flows, 
induced  effects,  tare  forces,  etc.),  accurate  drag  measurement  seemed  to 
be  the  most  difficult  problem  to  contend  with  in  designing  a  force 
balance.  The  solution  to  this  problem  seemed  to  lie  in  traversing  a 
pressure  probe  through  the  wake  of  a  given  airfoil  at  a  spanwise 
position  such  that  end  wall  and  tip  gap  effects  are  negligible. 
Integration  of  the  pressure  readings  would  then  give  an  accurate  section 
drag  coefficient.  Otsuka  and  Sugiyama's  work  seemed  to  indicate  that 


this  was  possible. 
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DRAG  COEFFICIENT,  CD 

Polar  diagram  for  the  NACA  4412  hydrofoil  in  nonca vitating 
flow.  The  data  points  are  from  the  present  tests.  The 
curves  are  from  the  following  sources: 

A  and  B.  Langley  Two-Dimensional,  Low  Turbulence 
Wind  Tunnel,  1949 

C.  J.  W.  Daily,  Hydrodynamics  Laboratory, 
C.I.T.,  1944. 

D  and  E.  NACA  Variable  Density  Wind  Tunnel,  1937. 

Figure  7.  Lif t-versus-drag  results  reported  by  Kermeen. 
Reproduced  from  Reference  (l)  . 


CHAPTER  III 


METHODS  AND  PROCEDURES 


3.1  Test  Program 

A  test  program  was  formulated  after  analyzing  all  the  available 
background  material.  The  first  step  would  be  to  develop  competence  with 
strain  gage  measurements  and  basic  wind  tunnel  techniques.  It  was  felt 
that  since  the  measurement  of  pitching  moment  coefficient  has  never 
proven  to  be  particularly  difficult  only  lift  and  drag  measurements 
would  be  considered.  This  would  also  simplify  the  hardware  manufactur¬ 
ing  problem.  An  existing,  doubly  supported  balance*  would  be  used  dur¬ 
ing  this  preliminary  phase.  This  balance  (Figure  8)  would  enable  forces 
to  be  measured  on  a  well-documented  airfoil  shape  in  the  absence  of  an 
incoming  boundary  layer. 

The  next  step  would  be  to  build  a  cantilever  balance  for  use  in  the 
48- in.  wind  tunnel.  This  balance  would  duplicate  as  closely  as  possible 
the  geometry  and  principle  of  operation  of  the  existing  12-in.  water 
tunnel  balance.  This  balance  would  be  used  to  measure  lift  coefficients 
for  the  selected  airfoil  profile  at  various  angles  of  attack.  For  these 
same  angles  of  attack,  a  pressure  traverse  would  be  used  to  measure  drag 


*A  doubly  supported  balance  refers  to  a  configuration  where  the  model  is 
supported  from  both  sides  of  the  tunnel  and  each  support  has  its  own 
force  measuring  device.  By  contrast,  a  cantilever  balance  herein 
refers  to  a  configuration  where  the  model  is  supported  from  only  one 
side  of  the  tunnel  with  a  single  force  measuring  device  attached  to  the 
support.  The  term  cantilever  balance  is  also  used  to  refer  to  the 
force  balance  described  in  Section  3.4. 
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DUMMY 


Figure  8.  Doubly  supported  balance.  Balance  elements  are  on 
either  side  of  tunnel.  Measurements  are  sunned  to 
produce  total  forces  on  model. 
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coefficients.  The  resulting  and  values  would  be  combined  to  pro¬ 
duce  a  C  -versus-C,  curve.  This  curve  would  be  verified  by  the  NACA 
x  d 

data  to  provide  a  data  base  for  determining  correction  factors  in  the 
final  force  balance  phase  of  the  investigation. 

The  next  step  would  be  a  variation  of  parameters  to  determine  their 
effect  on  measurements.  Airfoils  of  different  aspect  ratio,  disks  of 
various  sizes,  and  boundary  layers  of  varying  thickness  would  be  stud¬ 
ied.  These  data  would  be  analyzed  and  final  conclusions  drawn  regarding 
what  factors  significantly  influence  this  type  of  testing. 

3 . 2  Doubly  Supported  Balance  Phase 

The  doubly  supported  balance  as  depicted  in  Figure  8  was  designed 
to  measure  forces  on  an  airfoil  in  the  absence  of  an  incoming  boundary 
layer.  The  dummy  portions  of  the  airfoil  extend  into  the  flow  beyond 
the  thickness  of  the  sidewall  boundary  layers.  The  three-dimensional 
flows  occurring  at  the  airfoil-tunnel  wall  intersection  act  upon  the 
dummy  portion  of  the  foil  ends,  and  these  spurious  reactions  are  not 
transmitted  to  the  balance. 

The  force  cubes  were  mounted  on  two  parallel  wooden  (birch  veneer 
particle  board)  4  ft  x  8  ft  panels.  The  wind  tunnel  test  section  had  an 
octagonal  shape  and  installation  of  the  parallel  panels  converted  the 
test  section  to  a  rectangular  channel.  The  airfoil  setup  is  shown  on 
Figure  9.  Here  the  system  is  mounted  on  a  box- type  frame  used  to  cali¬ 
brate  the  balance  outside  of  the  wind  tunnel.  The  parallel  walls  shown 
on  Figure  9  simulate  the  panels  which  were  placed  in  the  wind  tunnel. 

The  force  cubes  were  mounted  on  the  outside  of  the  panels  (see 
Figure  8) . 


Figure  9.  Doubly  supported  balance  setup  in  calibration  box. 
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There  are  four  force  cubes  in  the  system  -  two  on  each  side  of  the 
tunnel.  Each  cube  is  designed  to  deflect  in  only  one  of  the  three  car¬ 
dinal  directions  and  to  be  rigid  in  the  other  two  directions.  The  cubes 
are  mounted  on  each  side  of  the  channel  so  that  one  cube  is  free  to 
deflect  in  the  horizontal  or  drag  direction  and  the  other  cube  is  free 
to  deflect  in  the  vertical  or  lift  direction.  Strain  gages  are  attached 
to  the  deflecting  faces  of  the  cubes.  The  strain  gages  produce  a  volt¬ 
age  which  is  proportional  to  the  force  applied  to  the  cube. 

Reference  (7)  provides  details  on  the  operation  of  strain  gages. 

The  cubes  are  mounted  in  a  fixed  position  on  the  tunnel  walls  such 
that  the  deflecting  axis  of  one  pair  is  perpendicular  to  the  tunnel 
flow.  The  shaft  which  transmits  forces  from  the  airfoil  to  the  force 
cubes  is  fixed  to  the  airfoil  with  four  setscrews.  The  angle  of  attack 
can  be  changed  by  loosening  the  setscrews,  rotating  the  airfoil  about 
the  shaft  to  the  desired  angle,  and  then  tightening  the  setscrews. 

The  cubes  were  calibrated  by  first  mounting  them  on  a  box  as  on 
Figure  9.  Known  weights  were  suspended  from  the  force  transmitting 
shaft  and  the  corresponding  output  was  recorded.  A  calibration  constant 
in  lb /volt  was  then  determined  for  each  cube.  After  the  lift  constants 
were  determined,  the  box  was  rotated  90°  so  that  the  drag  constants 
could  be  determined.  Additionally,  interactive  constants  were  deter¬ 
mined  by  loading  various  grams  of  "lift"  on  the  shaft  and  reading  the 
corresponding  volts  of  "drag.”  This  was  done  for  each  pair  of  cubes. 

An  interactive  constant  for  volts  of  lift  per  gram  of  drag  was  not 
determined  because  the  high  ratio  of  lift  to  drag  for  airfoils  (and  the 
design  of  the  force  gages)  makes  the  drag  effect  on  lift  negligible. 
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This  fact  was  verified  by  loading  the  cubes  with  representative  values 
of  lift  and  drag  simultaneously. 

Prior  to  mounting  the  balance  in  the  wind  tunnel,  a  pressure  probe 
was  traversed  vertically  in  the  channel  formed  by  the  two  wooden  panels 
to  verify  flow  uniformity.  These  results  are  shown  in  Figure  10. 
Traverses  were  also  made  perpendicular  to  the  sidewall  of  the  tunnel  and 
the  tunnel  ceiling  to  determine  boundary  layer  thickness.  These  results 
are  shown  in  Figures  11A  and  11B.  These  runs  were  all  made  at 
Re  =  330,000.  This  Re  was  chosen  since  it  was  the  highest  achievable  Re 
for  which  data  were  published  in  Reference  (3) .  The  plots  show  that 
there  was  a  uniform  flow  in  the  channel  and  that  the  sidewall  boundary 
layer  thickness  was  less  than  the  span  of  the  dummy  airfoils. 

Lift  and  drag  data  were  measured  using  the  doubly  supported  balance 
and  a  NACA  0012  airfoil.  The  airfoil  was  made  of  mahogany  and  had  a 
chord  of  6  in.  and  span  of  14.725  in.  Final  corrected  data  are  shown  in 
Figures  12  and  13.  The  corrections  applied  were  for  solid  blockage, 
wake  blockage,  lift  effect,  and  horizontal  buoyancy.  An  additional  cor¬ 
rection  was  applied  to  account  for  the  air  being  drawn  through  the  slot 
between  the  dummy  section  and  the  tested  airfoil.  In  this  small  gap, 
the  bar  which  transmits  forces  to  the  balance  is  exposed.  As  the  air¬ 
foil  becomes  more  heavily  loaded  at  higher  angles  of  attack,  the  larger 
pressure  difference  between  the  suction  and  pressure  surfaces  of  the 
airfoil  causes  a  higher  dynamic  pressure  to  be  sensed  by  the  bar.  A 
correction  term  for  this  additional  drag  was  determined  and  applied  to 
the  data.  This  correction  procedure  is  discussed  in  Section  3.3. 

To  determine  the  accuracy  of  the  balance  readings,  drag  was  meas¬ 
ured  by  traversing  a  total  pressure  probe  through  the  airfoil  wake  and 
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Figure  11B.  Velocity  profiles  measured  near  tunnel  ceiling 
The  traversable  probe  was  mounted  in  the  front 
and  rear  of  the  channel  as  shown  to  indicate 
boundary  layer  growth. 


-o-  NACA  ref.  data-,  Re  =  330,000 


—r*t*z  r  ~ 


rr 


33 


integrating  the  results  as  discussed  in  Appendix  B.  A  traversible  total 
pressure  probe  was  mounted  in  the  tunnel  one  chord  length  behind  the 
airfoil.  A  reference  pitot-static  tube  was  likewise  mounted  to  a  tunnel 
sidewall  at  the  same  downstream  location.  This  enabled  a  ratio  of  local 
to  free-stream  dynamic  pressures  to  be  obtained.  This  spatial  variation 
of  pressure  ratios  is  a  function  of  the  momentum  deficit  created  by  the 
model  and,  thus,  a  measure  of  the  section  drag  coefficient  at  the  par¬ 
ticular  spanwise  location  of  the  probe  traverse.  A  typical  wake  profile 
is  shown  in  Figure  14.  Wake  traverses  were  made  for  angles  of  attack  of 
approximately  0°,  4°,  and  6°.  The  spanwise  location  of  the  surveys  was 
also  varied.  Shown  in  Figure  15  are  the  results  of  the  various  surveys. 

For  direct  force  measurements  with  the  doubly  supported  balance, 
the  use  of  the  previously  mentioned  interactive  constants  proved  to  be 
erroneous.  The  original  mounting  technique  involved  mounting  the  force 
cubes  on  each  panel  using  a  carpenter's  level  to  assure  proper  vertical 
alignment.  Known  weights  were  then  placed  on  the  force  transmitting 
shaft  (simulating  a  "negative"  lift)  and  drag  voltages  were  noted.  This 
spurious  drag  voltage  would  then  be  subtracted  out  of  measured  drag 
voltages  during  each  test  as  a  step  in  the  data  reduction  process.  The 
interactive  constant  would  have  to  be  recalculated  each  time  the  balance 
was  remounted  or  adjusted  on  the  tunnel  wall.  The  results  of  this  type 
of  operation  were  unsatisfactory.  Lift  measurements  were  erratic  and 
drag  measurements  were  too  large  by  one  order  of  magnitude.  It  was 
ultimately  decided  to  eliminate  the  interactive  term  by  a  different 
mounting  technique.  The  cubes  were  adjusted  (with  the  voltmeter  oper¬ 
ating)  until  an  application  of  pure  negative  lift  (by  suspending  known 
weights  from  the  force  transmitting  shaft)  resulted  in  a  zero  (or 
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negligibly  small)  reading  on  the  drag  channel  of  the  voltmeter.  This 
technique  resulted  in  a  dramatic  improvement  in  the  accuracy  of  the 
read ings . 

Another  technique  that  required  attention  was  the  method  of  measur¬ 
ing  angle  of  attack.  It  became  obvious  that  the  angle  of  attack  had  to 
be  measured  as  precisely  as  the  various  force  coefficients.  The  ulti¬ 
mate  technique  for  measuring  the  angle  is  depicted  in  Figure  16.  A 
field  artillery  gunner's  quadrant  was  used  in  conjunction  with  an  air¬ 
foil  template.  The  template  was  constructed  so  that  it  contacted  the 
airfoil  at  the  leading  edge  and  at  two  other  points  on  the  airfoil's 
upper  surface.  The  upper  surface  of  the  template  was  machined  so  that 
when  the  device  was  placed  on  the  airfoil  the  machined  surface  was 
parallel  to  the  chordline.  The  gunner's  quadrant  measures  angles  from 
the  horizontal  (using  a  bubble  level)  in  mils  (6,400  mils  equals  360°). 
This  technique  enables  the  angle  of  attack  to  be  measured  accurately  to 
+0.056°. 

3.3  Experimental  Results  from  the  Doubly  Supported  Balance 

Shown  in  Figure  15  are  several  interesting  results.  For  a  given 
angle  of  attack,  the  section  drag  coefficient,  C^,  as  measured  by  the 
wake  traverse  increased  in  the  vicinity  of  the  slot  between  the  model 
and  the  dummy  airfoil  sections.  The  magnitude  of  the  increased  also 
appears  to  be  related  to  the  amount  of  lift  being  created  by  the  model. 
The  increased  drag  near  the  tip  of  the  airfoil  (slot  region)  is  quite 
similar  to  the  behavior  reported  by  Otsuka  and  Sugiyama  in 
Reference  (6) .  The  data  of  Figure  15  show  that  for  values  of  y/b 
greater  than  0.15  the  effect  of  the  tip  gaps  on  is  negligible.  As 
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Figure  16.  Angle-of-attack  measuring  technique  showing  artillery 
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the  plot  shows,  for  a  =  5.74°  the  section  measured  at  the  mid-span 
position  (C^  =  0.0119)  is  in  excellent  agreement  with  the  NACA  data 
(C^  =  0.0120).  The  area  under  the  curve  for  a  =  5.74°  was  measured  with 
a  planimeter  and  normalized  by  the  airfoil  span.  This  gave  a  computed 
=  0.0148,  which  is  in  good  agreement  with  the  balance  measured 
Cd  =  0.0140.  This  seems  to  indicate  that  the  balance  is  measuring  a 
total  drag  coefficient,  C^,  for  the  airfoil  which  includes  the  effect  of 
three-dimensional  flow  at  the  wing  tips.  The  wing  is  considered  to  be 
that  portion  of  the  two-dimensional  airfoil  model  excluding  the  two 
dummy  airfoil  segments. 

Another  interesting  result  is  the  shape  of  the  wake  behind  the  air¬ 
foil.  Figure  17  shows  how  the  wake  profile  takes  on  a  double  spiked 
shape  approximately  0.5  in.  inboard  of  the  airfoil  tip.  This  shape  is 
most  likely  due  to  shed  vorticity  and  further  indicates  an  absence  of 
two-dimensional  flow  conditions  in  this  vicinity. 

Lakshminarayana  (8)  conducted  studies  on  the  induced  drag  associ¬ 
ated  with  a  chordwise  slot  at  the  center  span  position  of  a  single  air¬ 
foil.  He  suggests  the  additional  drag  coefficient  to  be  given  by 
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where  =  lift  coefficient,  AR  =  aspect  ratio,  x  *  width  of  slot, 
c  »  airfoil  chord,  and  =  increment  of  additional  drag  coefficient 
due  to  slot  in  airfoil. 

A  further  correction  was  derived  to  account  for  the  flow  through 
the  slot  physically  Impinging  on  the  force  transmitting  shaft.  Assuming 
that  flow  through  the  slot  is  directed  at  a  45°  angle  to  the  free  stream 
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Figure  17.  Wake  profiles  in  vicinity  of  slot  for  doubly  supported 
balance  showing  development  of  double-spiked  profile. 
Also  shown  are  approximate  spanwise  locations  of 
traverses. 
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(flow  visualization  studies  conducted  by  Lakshminarayana  (8)  indicate 
that  this  is  a  reasonable  assumption) ,  an  additional  component  of  drag 
due  to  dynamic  pressure  on  the  bar  would  be  sensed  by  the  balance. 

ACp  =  cos  6(nondimensional  dynamic  pressure  on  bar)  ,  (2) 

where  0  is  the  angle  of  the  resultant  force  (Figure  18) .  Setting  for 

a  cylinder  (bar)  equal  to  1.0,  the  drag  component  of  force  exerted  on 
the  bar  becomes 

1  2 

A(Drag)^ar  =  cos  0 (-^  pV  )  (circumference  of  bar)  (gap  width)  (3) 

or 

A(Drag)bar  =  2  cos  6(-|-  pV2)(ird)t  ,  (4) 

where  d  =  diameter  of  shaft  and  t  =  width  of  gap.  The  factor  2  is 
included  because  there  are  two  slots,  one  on  each  end  of  the  model.  The 
work  of  Gearhart  (9)  suggests  that  flow  through  a  narrow,  sharp  edged 
slot  such  as  this  should  be  modified  to  reflect  an  effective  flow 
velocity  rather  than  a  theoretical  velocity.  Equation  (4)  then  becomes 

A(Drag)bar  =  cos  0p (0.8V)2iTdx  ,  (5) 

where  V  =  the  velocity  flowing  through  the  gap.  There  are  several  fac¬ 
tors  affecting  the  magnitude  of  this  velocity:  the  pressure  difference 
between  the  upper  and  lower  openings  of  the  slot  due  to  the  airfoil's 
suction  and  pressure  sides,  the  incoming  free-stream  velocity,  and  vis¬ 
cous  effects  that  occur  with  flow  through  a  narrow  slot.  After 
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consideration  of  these  effects,  the  magnitude  of  the  velocity  in  the  gap 
was  assumed  to  be  equal  to  the  free-stream  velocity.  This  value  was  a 
constant  for  all  airfoil  angles  of  attack  and  therefore  the  magnitude  of 
the  force  vector  acting  on  the  bar  was  likewise.  As  shown  in  Figure  18, 
when  a  increases,  the  angle  0  that  the  force  vector  makes  with  the  free- 
stream  decreases  and  the  drag  component  of  the  force  vector  increases. 

In  this  manner,  a  drag  component  is  created  which  varies  with  a  and 
free-stream  velocity. 

Equations  (1)  and  (5)  were  used  to  calculate  a  total  added  drag  due 
to  the  chordwise  slot  and  the  bar.  This  added  drag  was  subtracted  from 
the  balance  readings  for  drag  and  corrected  coefficients  are  shown  on 
Figure  13.  Also  shown  on  the  figure  are  three  values  as  measured  by 
the  wake  traverse  method.  The  traverses  were  done  at  the  mid-span  loca¬ 
tion  of  the  airfoil  for  three  different  a's.  The  data  of  Figure  13 
indicate  that  wake  traverses  for  this  size  airfoil  can  give  accurate  Cd 
measurements.  Further,  this  particular  balance  can  be  used  to  make 
accurate  two-dimensional  force  measurements  by  applying  suitable  cor¬ 
rection  factors. 

A  final  result  of  this  phase  of  the  program  was  the  realization  of 
the  importance  of  model  construction.  The  airfoil  was  examined  and  was 
found  to  have  a  varying  thickness  along  the  span.  The  thickness  tapered 
0.020  in.  from  tip  to  tip.  Figure  15  reflects  this  fact  in  that  the 
section  drag  coefficient  gently  decreases  from  mid-span  to  the  point 
near  the  tip  where  slot  effects  begin  to  take  over.  There  was  also 
found  to  be  a  1°  twist  in  the  model  from  tip  to  tip.  Since  all  angle  of 
attack  measurements  were  made  from  the  "high"  side  of  the  model,  0.5 
was  subtracted  from  each  measured  a  to  give  an  average  total  airfoil  a. 
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3 . A  Cantilever  Balance  Phase 

The  next  phase  of  the  program  concentrated  on  duplicating  the  con¬ 
ditions  encountered  in  the  water  tunnel.  A  two-component  cantilever 
balance  was  designed  and  built  which  was  similar  in  principle  to  the 
water  tunnel  balance.  The  balance  and  airfoils  are  shown  in  Figures  19, 
20,  and  21.  The  balance  used  the  principle  of  a  tension  member  as 
devised  by  Gurney  (lo) .  The  tension  member  is  a  beryllium-copper  strip 
of  approximately  0.006  in.  thickness.  The  functioning  of  the  tension 
strip  is  illustrated  in  Figure  22.  The  outer  two  clamps  fix  the  ribbon 
to  one  cube.  The  ribbon  was  clamped  while  under  a  tensile  force  greater 
than  one-half  the  largest  force  that  it  would  have  to  sense  (in  this 
case,  30  lb).  The  center  clamp  is  fixed  to  the  other  force  cube  and  is 
free  to  deflect  in  the  direction  of  the  long  axis  of  the  ribbon.  As  the 
cube  which  is  fixed  to  the  center  clamp  deflects  under  an  applied  force, 
one-half  of  the  tension  member  experiences  a  greater  tension  and  the 
other  half  experiences  less  tension  because  of  the  preload.  The  strain 
gages  sense  this  and  give  a  corresponding  voltage  which  is  calibrated 
for  force.  A  second  ribbon  is  clamped  to  a  second  cube  at  a  90°  angle 
to  the  first  ribbon.  The  two  force  cubes  effectively  have  one  face  in 
common.  This  design  enables  two  forces  to  be  measured  which  are  acting 
normal  to  each  other  (i.e.,  lift  and  drag).  The  advantage  of  this  sys¬ 
tem  is  that  the  cubes  require  very  little  clearance  space  to  accommodate 
the  small  deflections.  Tension  members  can  also  withstand  much  greater 
forces  than  the  strain  gage  cube  setup  (of  the  doubly  supported  balance 
system)  before  deflection  of  balance  parts  becomes  a  problem.  The 
resulting  high  sensitivity  balance  can  then  be  designed  very  compactly. 
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Figure  19.  Cantilever  balance  installed  on  tunnel  wall. 

Also  shown  are  angle-of-attack  indicator 
dial  and  wooden  square  frame  which  supports 
a  cover.  The  cover  shields  the  balance  from 
air  flowing  between  the  two-dimensional 
channel  wall  and  the  octagonal  wind  tunnel 
wall . 
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Figure  20.  Airfoil  with  18-in.  chord  and  no  disk  mounted  in  wind  tunnel. 


Figure  21.  Airfoil  with  9-in.  chord  and  11-in. -diameter  disk. 
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The  balance  and  airfoil  system  is  shown  in  simplified  form  in 
Figure  23.  Originally,  the  balance  was  mounted  so  that  the  ribbons 
formed  a  45°  angle  with  the  free  stream.  This  configuration  was  tried 
in  an  attempt  to  have  each  ribbon  measure  approximately  the  same  gross 
force.  In  the  doubly  supported  phase,  the  drag  cubes  measured  forces  of 
a  lower  order  of  magnitude  than  the  lift  cubes.  It  was  theorized  that 
this  difference  in  the  order  of  magnitude  of  the  two  forces  could  have 
bein  a  source  of  the  problem  with  the  water  tunnel  balance.  Ultimately, 
for  reasons  detailed  in  Section  3.5,  the  balance  was  aligned  so  that  one 
ribbon  measured  normal  forces  and  the  other  ribbon  measured  axial  forces 
as  depicted  in  Figure  23.  The  balance  rotated  with  the  model  as  angle 
of  attack  was  changed. 

Two  airfoil  models  of  different  AR  were  then  fashioned  for  use  with 
the  cantilever  balance.  A  NACA  0012  airfoil  was  machined  of  aluminum. 
This  model  has  a  chord  of  9  in.  and  a  span  of  18.375  in.  for  an  aspect 
ratio  of  2.04.  The  model  weighed  approximately  11  lb  14  oz.  A  second 
model  was  made  of  a  two-part  expandable  urethane  foam.  A  master  airfoil 
was  machined  from  aluminum  (span  =  18.375  in.,  chord  =  18  in., 
weight  =  40  lb).  A  high-strength  mold  was  then  made  from  this  master. 
The  two-part  foam  was  mixed  and  poured  into  the  mold.  The  foam  expanded 
within  the  mold  and  gave  a  lightweight  duplicate  airfoil.  Some  shrink¬ 
age  resulted  from  this  manufacturing  technique  (approximately  0.001  in. 
per  in.),  but  the  finished  product  was  considered  quite  satisfactory. 

The  final  model  weight  was  approximately  10  lb  12  oz. 

The  calibration  procedure  for  the  cantilever  balance  was  similar  to 
that  of  the  doubly  supported  balance.  Each  ribbon  was  calibrated  with 
known  weights  after  the  balance  was  mounted  in  the  tunnel.  The 


Figure  23.  Simplified  sketch  of  cantilevered  balance.  Balance 
mounts  on  tunnel  wall.  A  shaft  protrudes  through 
the  wall  and  is  fixed  to  the  tested  airfoil.  The 
balance  rotates  with  the  airfoil. 
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calibrations  were  performed  after  the  interactive  constants  were  elec¬ 
tronically  zeroed  out  in  a  fashion  similar  to  that  of  the  previous 
phase.  Reduction  of  the  raw  data  for  this  balance  configuration  is 
explained  in  Appendix  C. 

The  first  survey  made  in  this  test  phase  was  a  measure  of  the  tun¬ 
nel's  static  pressure  gradient  to  determine  the  horizontal  buoyancy  cor¬ 
rection.  Static  pressure  taps  were  placed  in  one  of  the  4  ft  x  8  ft 
panel  walls  at  four  axial  locations.  The  taps  were  located  1,  3,  5,  and 
7  ft  from  the  panel's  leading  edge  and  1  ft  from  the  panel  floor.  The 
results  are  shown  in  Figure  24.  The  airfoil  was  not  mounted  in  the 
tunnel  for  these  tests.  Pressures  are  plotted  versus  channel  position 
for  four  velocities  as  indicated.  The  pressures  were  nondimensionalized 
by  dynamic  pressure  to  give  an  average  value  of  (dP/dJl)/q  =  0.01236  ft  \ 
where  dP/dS.  =  static  pressure  change  per  unit  length  in  channel 

(lb/ft2)/ft  and  q  =  (l/2)pV  2. 

o 

The  balance  was  mounted  28  in.  from  the  leading  edge  of  the  wooden 
panels.  This  position  was  chosen  to  ensure  a  thin  incoming  boundary 
layer.  It  was  felt  that  the  boundary  layer  could  be  thickened  later  for 
comparison  testing.  However,  adding  2  in.  of  1/100  sandpaper  to  the 
leading  edges  of  the  panels  failed  to  make  a  significant  change  in  the 
boundary  layer  thickness.  After  further  consideration,  it  was  felt  that 
the  balance  should  have  been  mounted  farther  downstream  on  the  channel 
wall.  This  would  have  given  leading  edge  disturbances  more  time  and 
distance  to  grow  into  a  thicker  boundary  layer.  Results  of  the  boundary 
layer  surveys  are  shown  on  Figure  25. 

The  next  step  in  the  program  was  to  mount  the  airfoil  with  the 


chord  of  9  in.  An  11-in. -diameter  disk  was  also  mounted  on  one  end  of 


total 
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the  airfoil  simulating  the  spindle  disk  configuration  of  Kermeen.  With 
the  airfoil  in  place,  wake  measurements  were  made  with  the  traversable 
probe  as  in  the  two-component  balance  phase.  Several  of  the  wake  pro¬ 
files  are  shown  in  Figures  26  and  27.  The  profiles  show  the  double 
spike  in  the  vicinity  of  the  airfoil-wall  intersection  similar  to  that 
seen  with  the  doubly  supported  balance.  Figure  28  contains  a  plot  of 
section  drag  coefficients,  C^,  versus  span  location.  The  plots  also 
indicate  that  beyond  a  certain  spanwise  location  (y/b  >  0.22  in  this 
case)  the  effect  of  the  tunnel  walls  on  the  probe  Cd  measurements  is 
negligible. 

The  9-in.  airfoil-ll-in.  disk  combination  was  then  tested  at  vari¬ 
ous  angles  of  attack.  Next  the  disk  by  itself  was  mounted  to  the  bal¬ 
ance.  The  airfoil  was  mounted  from  the  opposite  side  of  the  tunnel  with 
a  gap  between  the  airfoil  tip  and  disk  of  0.004  in.  In  this  manner, 
tare  forces  on  the  disk  were  determined.  The  tare  forces  were  then 
subtracted  from  the  total  measured  forces.  The  total  forces,  tare 
forces,  and  corrected  data  are  shown  in  Figures  29,  30,  and  31.  The 
corrected  data  include  the  traditional  corrections  previously  mentioned 
in  connection  with  Kermeen' s  work. 

The  lift  measurements  obtained  with  the  balance  were  combined  with 
drag  data  measured  with  the  traversable  probe  (for  given  angles  of 
attack)  and  plotted  in  comparison  with  the  NACA  data,  Figure  31.  There 
is  excellent  agreement  between  the  two  sets  of  data.  These  two  sets  of 
data  established  a  reference  base  for  the  determination  of  correction 
factors. 

At  this  point,  there  existed  an  error  in  the  balance-measured 
which  varied  from  zero  to  approximately  30%  of  the  NACA  reference  data 
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Figure  26.  Wake  profiles  in  the  vicinity  of  the  airfoil-disk 
intersection  showing  growth  of  a  second  wake 
behind  airfoil  upper  surface. 


y / B  (2 B  =  tunnel  width) 


momentum  deficit  technique. 
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at  the  higher  values  (Figure  29).  The  error  trend  was  similar  in 
nature  to  that  reported  by  Daily  and  Kermeen.  The  first  attempt  to 
deduce  the  cause  of  this  error  involved  studying  the  effect  of  the  gap 
between  the  airfoil  tip  and  the  channel  wall. 

A  fiber  optics  instrument  (Fotonic  Sensor)  was  used  to  measure  pre¬ 
cisely  the  size  of  the  airfoil  tip  gap.  A  probe  capable  of  directing  or 
detecting  a  small  beam  of  light  was  imbedded  in  the  tunnel  wall.  The 
light  would  be  reflected  off  the  aluminum  airfoil  tip  and  give  an  accu¬ 
rate  measure  of  the  distance  between  the  probe  (flush  with  tunnel  wall) 
and  airfoil  tip.  Using  this  technique,  the  gap  was  varied  from 
0.001  in.  to  0.010  in.  while  maintaining  an  angle  of  attack  of  5.9°. 

For  this  range  of  end  gap,  it  was  found  that  there  are  no  significant 
effects  on  measurement  and  that  CL  values  decreased  by  less  than  1%. 
These  data  reinforce  previous  investigations  (ll)  that,  if  the  tip  gap 
can  be  held  sufficiently  small,  viscous  forces  predominate  and  the 
effect  of  the  gap  on  airfoil  forces  becomes  negligible. 

In  an  attempt  to  determine  the  importance  of  the  disk  size,  the 
9-in.  chord  model  was  tested  without  an  attached  disk.  The  results  are 
shown  in  Figures  32A  and  32B.  Data  taken  with  a  disk  attached  to  the 
9-in.  airfoil  and  tare  forces  subtracted  are  also  shown  for  comparison. 

tare  forces  due  to  the  disk  were  negligible  and  are  not  included  in 
Figure  32A.  The  traditional  buoyancy  and  blockage  corrections  have  been 
applied.  As  Figure  32B  shows,  the  disk  does  not  appear  to  be  essential. 
This  is  important  Information  from  a  manufacturing  standpoint. 

At  one  point,  the  C^-versus-a  curve  showed  a  peculiar  behavior. 

Near  C  ,  the  C  values  showed  evidence  of  a  sharp  stall  which  was  not 
max 

in  agreement  with  the  NACA  data  (the  open  squares  in  Figure  33) .  The 
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Figure  32A.  versus  a  (corrected)  for  18-in. 

9-in. -chord  airfoils  without  disks 
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Figure  33.  versus  a  for  9-in. -chord  aluminum  airfoil 

with  and  without  sandpaper  on  leading  edge. 
The  graph  indicates  that  laminar  separation 
causes  sharp  stall  characteristics. 
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smooth  finish  on  the  aluminum  model  suggested  that  possibly  laminar  sep¬ 
aration  was  occurring.  To  investigate  this,  a  1/4-in.  strip  of  #100 
sandpaper  was  glued  to  the  leading  edge  of  the  airfoil  to  trip  the 
boundary  layer.  The  C^-versus-a  curve  was  then  "smoothed  out"  to  the 
appropriate  shape  as  shown  by  the  solid  circles  in  Figure  33.  The  C 
values  were  lower  due  to  the  sandpaper,  but  the  behavior  of  the  curve  in 
the  stall  region  proved  that  laminar  separation  had  occurred. 

In  a  further  attempt  to  discover  the  source  of  additional  drag 
readings,  the  flow  on  the  spindle  side  of  the  airfoil  was  considered. 

It  was  felt  that  there  was  a  possibility  that  at  high  C^'s  the  flow 
could  be  drawn  through  the  gaps  between  the  model  tips  and  the  tunnel 
walls.  On  the  balance  side,  this  flow  could  be  directly  impinging  on 
the  force  transmitting  bar  and  causing  the  balance  to  read  slightly 
higher  values.  This  is  the  same  condition  considered  in  the  doubly 
supported  balance  phase.  To  find  out  if  this  force  was  significant,  a 
stub  spindle  arrangement  was  fashioned  as  depicted  in  Figure  34. 

The  test  setup  was  designed  to  measure  only  spurious  drag  readings. 
The  test  was  further  designed  to  determine  whether  this  type  of  addi¬ 
tional  drag  was  of  a  sufficiently  large  magnitude  to  explain  the  behav¬ 
ior  of  the  original  water  tunnel  balance.  The  airfoil  was  mounted  from 
the  opposite  wall  of  the  tunnel  as  in  the  previous  tare  tests.  A  stub 
spindle  which  protruded  0.020  in.  into  the  airstream  was  attached  to  the 
balance  (Figure  34).  A  continuity  tester  with  leads  attached  to  the 
airfoil  and  balance  was  used  to  keep  the  gap  between  the  stub  and  the 
airfoil  as  small  as  possible.  Under  these  conditions,  air  flows  in  the 
vicinity  of  the  model-wall  intersection  are  closely  duplicated  and  the 
balance  measured  only  the  forces  exerted  on  the  shaft.  The  drag  forces 
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Figure  34. 


Stub  spindle  setup  to  determine  forces  acting  on  force 
transmitting  shaft  by  flows  through  end  gap. 
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were  nondimensionalized  by  the  airfoil  planform  area  so  that  the  result¬ 
ing  increment  could  be  subtracted  from  the  model  C^.  The  results  are 
shown  in  Table  1. 

As  shown  in  Table  1,  small  amounts  of  drag  were  measured  by  the 
balance.  The  value  of  this  drag  increment  increased  with  angle  of 
attack.  Based  on  this  information,  it  was  then  felt  that  the  shaft- 
airfoil  junction  should  be  made  to  duplicate  more  closely  that  used  in 
the  water  tunnel.  The  hydrofoil-shaft  arrangement  used  in  the  water 
tunnel  is  shown  in  Figure  35.  At  the  point  of  attachment,  there  was  a 
0.090-in.  gap  between  the  plane  of  the  smooth  face  of  the  shaft  and  the 
hydrofoil  tip.  This  meant  that  the  force  transmitting  shaft  was  effec¬ 
tively  a  long  rectangular  shaped  beam.  It  was  theorized  that  this  broad 
beam  lying  in  the  end  gap  flow  could  result  in  large  measured  drag 
increments.  The  wind  tunnel  model  was  modified  as  shown  in  Figure  36. 

A  wooden  spacer  made  the  shaft  diameter  larger,  and  the  area  between  the 
larger  shaft  and  airfoil  was  waxed-in  to  duplicate  the  0.090-in. 
standoff  of  the  water  tunnel  model.  This  configuration  (Figure  36)  was 
tested  at  a  =  7°.  The  amount  of  spurious  drag  sensed  obviously  did  not 
cause  the  drag  measurement  to  vary  by  an  order  of  magnitude.  Therefore, 
the  flow  impinging  on  th  Horce  transmitting  shaft  was  ruled  out  as  a 
major  source  of  error. 

The  final  portion  of  the  experimental  phase  was  the  testing  of  the 
large  airfoil  (chord  -  18  in.).  This  test  was  conducted  without  a  disk 
on  the  balance  end  of  the  model.  The  results  are  shown  in  Figures  32A 


and  32B. 


TABLE  1 


Drag  Increment  Measured  by  Cantilever  Balance  on 
Stub  Spindle  of  Figure  34. 

Re  =  18,330  (based  on  spindle  dia  =  0.5  in.) 


a  (deg) 

C,  (based  on  c  -  9  in.) 
a 

0.0 

0.00001 

7.1 

0.0002 

7.77 

0.00055 

8.9 

0.0008 

10.8 


0.0013 
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Figure  35.  Water  tunnel  shaft  and  hydrofoil  as  used  in 
experiments.  Inset  illustrates  gap  that 
existed  on  balance  side  of  hydrofoil. 


WOOD  SPACER  TO 
EXTEND  THICK 


Figure  36.  Modification  applied  to  wind  tunnel  airfoil  to 
duplicate  configuration  of  water  tunnel 
hydrofoil  and  shaft  of  Figure  35. 
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3.5  Experimental  Results  for  Cantilever  Balance  Phase 

Results  from  this  phase  of  the  experiments  indicate  that  the  mount¬ 
ing  of  this  type  of  balance  is  critical.  At  the  higher  angles  of 
attack,  the  ratio  of  lift  force  to  drag  force  becomes  very  high.  The 
drag  force  appears  to  be  too  small  to  measure  as  the  difference  or  sum 
of  two  small  numbers.  When  the  balance  was  mounted  with  the  ribbon  at  a 
45°  angle  with  respect  to  the  free  stream,  the  axial  force  was  shared 
between  both  tension  members.  When  the  balance  was  remounted  so  that 
one  tension  member  measured  all  of  the  axial  force  and  the  other  tension 
member  measured  all  of  the  normal  force,  the  results  were  dramatically 
better . 

Another  significant  result  stems  from  the  Fotonic  Sensor  tests. 

The  data  indicate  that  within  a  certain  range  the  size  of  the  gap 
between  the  airfoil  tip  and  tunnel  wall  has  a  negligible  effect  on 

and  C,  measurements. 

a 

The  idea  of  using  a  disk  on  the  balance  end  of  the  model  yielded 
interesting  results.  Experimental  data  on  tests  with  and  without  the 
disk  indicate  that  the  inclusion  of  the  disk  is  not  necessary.  This  is 
from  an  aerodynamic  standpoint,  however.  Use  of  a  disk  to  facilitate 
making  angle-of-attack  changes,  aiding  in  alignment,  etc.,  may  be  desir¬ 
able  from  a  hardware  design  standpoint. 

Another  interesting  result  concerns  the  effect  of  aspect  ratio. 

The  airfoil  of  9-in.  chord  had  an  aspect  ratio  of  2.  The  18-in.  airfoil 
had  an  aspect  ratio  of  1.0.  From  the  experimental  data  (Figures  32A  and 
32B) ,  the  effect  of  aspect  ratio  variation  appears  to  be  small  (at  least 
for  these  two  AR's). 
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CHAPTER  IV 

ANALYSIS  AND  DISCUSSION 

As  the  experimental  results  demonstrated,  measuring  and  C^  sep¬ 
arately  (by  the  most  accurate  technique  for  each)  and  combining  the 
results  produced  a  satisfactory  reference  curve.  Data  generated  in  this 
manner  were  in  good  agreement  with  the  NACA  data.  In  order  to  make  the 
balance  system  an  accurate,  two-dimensional  force  measuring  device,  a 
correction  must  be  derived  to  account  for  the  difference  between  the 

reference  C  -versus-C,  curve  and  the  balance-produced  curve, 
x.  a 

The  difficulty  in  making  mechanical  balances  perform  accurately  is 
the  drag  measurement.  For  a  given  airfoil  operating  at  normal  angles  of 
attack,  drag  forces  tend  to  be  much  smaller  than  the  corresponding  lift 
forces.  In  both  the  water  tunnel  and  the  wind  tunnel,  when  the  balance 
was  set  up  to  measure  drag  or  axial  forces  by  combining  the  output  of 
two  measuring  elements,  the  drag  results  were  in  error  by  an  order  of 
magnitude.  Reorienting  the  wind  tunnel  balance  so  that  one  measuring 
element  measured  normal  forces  and  the  other  element  measured  axial  for¬ 
ces,  with  no  interactive  terms,  brought  the  drag  data  into  the  regime 
achieved  by  previous  authors.  A  simple  reorientation  of  the  balance  did 
not  achieve  the  same  results  for  the  water  tunnel  balance.  However,  the 
water  tunnel  balance  had  problems  with  visible  deflections  of  the  model 
when  the  balance  was  in  the  most  favorable  (from  a  drag  measurement 


standpoint)  orientation. 
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Assuming  the  proper  balance  orientation  had  been  selected,  there 
was  still  an  increment  of  drag  to  be  accounted  for.  An  increment  of 
drag  varying  from  zero  at  =  0  to  approximately  35%  (or  greater)  of 
the  published  NACA  value  at  higher  values  was  reported  by  Kermeen  and 
Daily  in  their  tests  of  a  NACA  4412  airfoil.  The  same  general  behavior 
of  was  found  in  the  results  of  the  present  report  for  a  NACA  0012 
airfoil.  Since  this  behavior  is  found  in  mechanical  balance  experiments 
only  and  mechanical  balances  measure  the  forces  created  by  non-two- 
dimensional  flows  on  the  model,  it  was  assumed  that  the  additional  C^ 
increment  documented  in  the  three  reports  was  due  to  non-two-dimensional 
flow  effects  on  the  model. 

It  has  long  been  known  that,  when  a  strut  intersects  a  flat  surface 
in  the  presence  of  a  nonuniform  flow,  secondary  flows  are  created.  In 
the  case  of  an  airfoil  completely  spanning  a  channel,  the  boundary  layer 
on  the  tunnel  sidewall  causes  a  so-called  horseshoe  vortex  to  be  formed 
which  engulfs  the  airfoil-wall  intersection.  See  Figure  37.  This 
horseshoe  vortex  is  a  region  of  "contaminated"  flow  in  that  two- 
dimensional  flow  conditions  do  not  exist  there.  In  order  to  correct  the 
balance  measurements  to  two-dimensional  forces  on  the  hydrofoil,  this 
contaminated  region  must  be  accounted  for. 

There  are  several  methods  of  dealing  with  secondary  flows  of  this 
nature.  Suction  or  blowing  devices  can  be  installed  which  will  physi¬ 
cally  remove  the  sidewall  boundary  layers.  Removing  the  boundary  layers 
eliminates  secondary  flows  and  enables  the  balance  to  measure  forces 
that  exist  on  a  model  immersed  in  purely  two-dimensional  flow.  Another 
approach  is  to  permit  the  model  to  experience  the  secondary  flows  and 
apply  a  theoretical  or  empirical  correction  to  the  data.  Since  suction 
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Figure  37.  Horseshoe  vortex  created  by  intersection  of  velocity 
gradient  and  curved  airfoil  surface  protruding  from 
tunnel  wall. 
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or  blowing  can  lead  to  excessive  cavitation  in  water  tunnels,  the  use  of 
a  correction  factor  is  preferred.  To  facilitate  deriving  a  correction 
factor,  some  insight  into  the  physical  mechanism  causing  the  additional 
drag  associated  with  secondary  flow  must  first  be  obtained. 

An  interesting  study  of  this  problem  was  recently  completed  by 
Barber  (l2) .  He  considered  the  additional  drag  created  by  a  strut  pro¬ 
truding  from  a  wall  for  various  incoming  boundary  layer  thicknesses.  He 
found  that  the  size  of  the  horseshoe  vortex  varied  directly  with  the 
thickness  of  the  incoming  boundary  layer.  He  also  found  that  the  por¬ 
tion  of  the  airfoil  where  flow  separation  occurred  varied  inversely  with 
the  size  of  the  horseshoe  vortex.  He  concluded  that  with  a  large  horse¬ 
shoe  vortex  viscosity  caused  high  momentum  fluid  to  be  entrained  in  the 
corner  where  the  airfoil  trailing  edge  and  wall  intersected  as  shown  in 
Figure  38.  This  influx  of  high  momentum  fluid  enables  the  flow  better 
to  withstand  the  adverse  pressure  gradient  existing  in  the  corner  and 
retards  flow  separation.  As  the  figure  shows,  a  thin  vortex  is  not  able 
to  entrain  as  much  of  the  high  momentum  fluid  and  a  larger  separated 
zone  exists.  Barber  made  numerous  flow  visualization  experiments  to 
verify  the  existence  of  this  flow  condition. 

Assuming  that  the  separated  region  near  the  trailing  edge  of  the 
airfoil-wall  intersection  is  responsible  for  the  increment  of  additional 
drag  measured,  one  can  link  a  correction  factor  to  the  thickness  of  the 
incoming  boundary  layer.  Support  for  this  idea  comes  from  the  work  of 
Hawthorne  (l3) . 

Hawthorne  derived  the  following  expression  for  the  energy  in 
* 

secondary  flows,  Dg  ,  created  by  strut-wall  intersections: 


PROPOSED  MODEL  OF  THICK  BOUNDARY-LAYER-STRUT  INTERACTION 


PROPOSED  MODEL  OF  THIN  BOUNDARY-LAYER-STRUT  INTERACTION 

Figure  38.  Barber's  model  for  the  flow  conditions  occurring 
in  the  vicinity  of  an  airfoil-tunnel  wall 
intersection.  Reproduced  from  Reference  (l2) . 
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and  p  =  density,  =  free-stream  velocity,  c  =  chord,  t  =  airfoil 
thickness,  and  n  =  4 (l  +  (1/2)  (t/c))  /  (l5tr (6/c) )  .  Hawthorne  shows  how 
f(n)  varies  with  boundary  layer  thickness  (6/c)  as  reproduced  in 
Figure  39.  The  plot  presents  data  from  a  bicusped  strut  profile  in  an 
exponential  boundary  layer  with  strut  thickness  values  of  0.05  and  0.25. 
The  figure  shows  that  f(n)  increases  with  6/c  to  a  maximum  value  at 
6/c  =  0.1.  Equation  (6)  states  that  the  energy  in  these  secondary  flows 
is  proportional  to  airfoil  thickness  to  the  fourth  power  and  reaches  a 
maximum  when  6/c  is  approximately  0.1.  Hawthorne's  analysis  was  based 
on  a  bicusped  airfoil  shape  and  an  exponential  boundary  layer.  Although 
the  theory  does  not  hold  for  all  airfoil  shapes  or  boundary  layer  pro¬ 
files,  it  is  probably  not  unfair  to  assume  that  generally  the  energy  in 
secondary  flows  for  this  type  of  airfoil-tunnel  wall  intersection  is 

De*  -  K(t/c)4(|^)  ,  (8) 


where  K  *  constant.  Assuming  further  that  the  ability  of  the  horseshoe 
vortex  to  entrain  high-momentum  flow  in  the  strut  intersection  and 
retard  separation  is  directly  related  to  the  energy  in  the  vortex  means 
that  the  drag  correction,  AC^,  becomes 
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f(C  ,  C  ,  a,  AR) 

AC  =  K  - - - 2 - —  .  (9) 

[(6/C)/0. l) (t/c)4 

The  C^,  a,  and  AR  effects  must  be  derived  from  an  analysis  of  the  pres¬ 
ent  experimental  data. 

Shown  in  Figures  AO-42  are  C^-versus-C^  data  as  reported  by  Daily, 

Kermeen,  and  the  present  study.  Also  shown  is  the  deviation  between  the 

experimental  data  and  the  published  NACA  data.  This  deviation  is  the 

required  correction  to  make  experimental  data  match  the  NACA  published 

data.  The  three  sets  of  data  have  several  similarities.  The  deviation 

of  drag  is  zero  at  =  0  for  all  three  investigations.  The  AC^  curves 

for  all  three  experiments  increase  to  the  point  where  dC^/da  is  no 

longer  constant  and  then  AC^  decreases.  If  (dC^/da)^  represents  the 

slope  of  the  lift  curve  in  the  constant  slope  portion  of  the  curve,  then 

the  AC,  curve  seems  to  behave  as 
o 


dC  /da 

1Cd  •  ^(d-Jdc.)  > 

K  O 


(10) 


for  all  three  experiments. 

Assuming  Equation  (10)  represents  the  contribution  of  a  to  the  drag 

correction,  one  sees  that  AC,  now  becomes 

d 
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Figure  40.  Daily's  results  (C^  versus  C^)  compared  with  NACA 

published  data  for  a  NACA  4412  airfoil.  Also 
shown  is  the  value  of  AC^  required  to  correct 

Daily's  curve  to  the  NACA  curve. 
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Figure  41.  Kermeen's  results  (C.  versus  C.)  compared  with  NACA 

a 

published  data  for  a  NACA  4412  airfoil.  Also  shown 
is  the  value  of  AC^  required  to  correct  Kermeen's 

curve  to  the  NACA  curve. 
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Figure  42.  ARL/PSU  results  (C.  versus  C.)  compared  with  NACA 

Jt  a 

published  data  for  a  NACA  0012  airfoil.  Also 
shown  is  the  value  of  required  to  correct  the 

ARL/PSU  curve  to  the  NACA  curve. 
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where  C£  =  dC£/da  and  C£  =  (dC  /da)  .  As  the  experimental  data  show, 

a  a  0 

o 

the  drag  correction  must  be  zero  at  C£  =  0.  Thus,  the  C£  contribution 

should  be  directly  proportional  to  the  AC,  term.  The  assumed  C,  cor- 

d  d 

rection  can  be  the  term  which  "individualizes"  the  correction  to  spe¬ 
cific  airfoils  if  it  is  set  equal  to  C,  (minimum  drag  coefficient)  of 

d 

o 

the  tested  airfoil.  The  correction  now  takes  the  form 


Cd  /C,  >1/2 

o  a  a 

AC  =  K  - - — r-  f  (AR) 

((S/c)/0.l) (t/c)4 


(12) 


At  this  point.  Equation  (12)  can  be  applied  to  the  data  of 
Figures  40-42.  On  a  trial-and-error  basis,  setting  K  *  0.00015  and 
f(AR)  =  1/i/AR  ,  one  finds  for  the  final  correction  that 


ACd 


0.00015(Cd  C£(C£  /C  )1/2) 
o  a  a 

_ o _ 

((6/c)  /O.l)  (t/c)4^R(l/2) 


(13) 


Shown  in  Figures  43-45  are  the  results  of  the  three  experiments 
with  Equation  (13)  applied  to  the  data.  The  correction  works  extremely 
well  for  the  ARL/PSU  data.  There  is  also  exceptionally  good  agreement 
when  it  is  applied  to  Kermeen's  data.  With  Kermeen's  data,  there  is  a 
slight  deviation  at  the  higher  C£  values.  It  must  be  noted  that  in  the 
range  C£  •  1.1  Kermeen's  curve  exceeds  the  published  NACA  data  by 
approximately  90%.  Another  consideration  is  the  advances  made  in 
instrumentation  between  Kermeen's  tests  and  the  present  study.  Also 
noteworthy  is  the  fact  that  the  Reynolds  numbers  of  the  two  curves 


Figure  44.  Kermeen's  data  corrected  by  Equation  (13) 
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a  18"  FOIL,  NO  DISK,  AR  =  1.0,  6/c  =  0.083 
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differ  by  88,000.  Higher  Reynolds  numbers  tend  to  make  the  versus 

curve  "flatten  out."  Daily's  data  show  excellent  agreement  even  in  the 

high  regime.  It  is  interesting  to  note  that  the  maximum  error  in 

Daily's  uncorrected  data  is  less  than  50%  of  the  NACA  data  or  almost 

half  of  the  error  reported  by  Kermeen. 

A  feature  of  Equation  (13)  which  is  of  critical  importance  is  the 

fact  that  all  the  terms  of  the  equation  are  easily  determined  with  the 

existing  balance  system.  It  has  been  demonstrated  that  C^,  ,  and 

a 

are  all  obtainable  with  the  present  balance  using  only  the  tradi- 
a 

o 

tional  boundary  corrections.  The  only  drag  dependent  term  in 

Equation  (13)  is  C,  .  C,  is  the  one  value  on  the  C  -versus-C  curve 
a  a  Jc  a 

o  o 

that  the  balance  can  measure  and  requires  no  correction.  Had  the  cor¬ 
rection  been  a  function  of  the  airfoil's  actual  or  corrected  drag  coef¬ 
ficient,  such  as  the  term,  it  would  be  of  no  practical  value  since 

a 

for  a  new  airfoil  design  the  corrected  values  are  unknowns. 


CHAPTER  V 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

5.1  Summary  and  Conclusions 

With  the  derivation  of  Equation  (13),  all  five  of  the  previously 
stated  objectives  have  been  met.  The  first  stated  objective  was  to 
determine  what  caused  the  original  water  tunnel  balance  to  fail  to  meas¬ 
ure  correctly  the  C^  values.  When  measurements  were  begun,  both  bal¬ 
ances  (water  and  wind  tunnel)  were  oriented  in  such  a  manner  that  they 
measured  drag  by  summing  or  differencing  two  small  force  measurements. 

In  both  cases,  the  resulting  C^  values  were  too  large  by  an  order  of 
magnitude.  When  the  wind  tunnel  balance  was  reoriented  such  that  one 
measuring  element  sensed  the  total  axial  force  on  the  model  and  the 
other  element  sensed  the  total  normal  force,  the  resulting  Cd  values 
were  within  the  regime  of  measurements  reported  by  previous  authors. 
Similar  reorientation  of  the  water  tunnel  balance  did  not  produce  the 
same  results.  However,  when  the  water  tunnel  balance  was  reoriented, 
the  flexures  between  the  measuring  elements  allowed  the  hydrofoil  to 
deflect  substantially  under  load,  rendering  the  reading  suspect.  As 
Figures  5A  and  5B  show,  both  C^  and  C^  data  had  a  large  error.  This 
proof  is  not  absolute,  but  in  light  of  the  failure  of  all  subsequent 
attempts  to  duplicate  various  error  causing  conditions  in  the  wind  tun¬ 
nel  the  above  reason  emerges  as  the  most  probable  explanation. 

The  second  objective  was  to  build  and  test  a  wind  tunnel  balance 
which  was  similar  in  principle  to  the  water  tunnel  balance.  This  was 
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accomplished  and  the  measured  results  were  as  good  as,  if  not  better 
than,  the  data  of  Kermeen  (l)  and  Daily  (5). 

The  third  objective  was  to  determine  what  factors  affect  flow  past 
a  hydrofoil  in  a  two-dimensional  channel  or,  just  as  important,  what 
factors  have  a  negligible  effect. 

The  experimental  results  show  that  several  factors  are  either  neg¬ 
ligible  or  correctable.  The  effect  of  a  small  gap  between  the  airfoil 
tip  and  channel  wall  is  negligible  provided  it  is  kept  within  a  certain 
range.  Gap  sizes  up  to  0.020  in.  for  the  18-in. -chord  model  show  no 
significant  effect.  Practically,  if  the  gap  can  be  kept  as  small  as 
possible  from  a  manufacturing  tolerance  standpoint  (gap/chord 
ratio  _<  0.002  in  this  study),  it  is  unlikely  that  any  significant  error 
will  be  introduced  by  the  existence  of  an  end  gap.  Another  factor 
deemed  noncritical  is  flow  in  the  vicinity  of  the  force  transmitting 
shaft.  Use  of  a  disk  on  the  shaft  end  of  the  airfoil  to  protect  the 
shaft  from  spurious  forces  does  not  appear  to  be  necessary.  The  varia¬ 
tion  of  aspect  ratios  used  in  the  current  test  program  was  limited 
(AR  =  1.0,  2.0,  and  3.3).  The  effect  of  varying  AR  between  1.0  and  2.0 
is  negligible  or  at  least  difficult  to  discern.  However,  when  AR  =  3.3 
is  considered,  as  in  Daily's  data,  there  is  an  effect  and  the  drag  cor¬ 
rection  factor  needs  modification  to  reflect  this.  Experimentally, 
f(AR)  ■  1//AR  serves  as  a  suitable  approximation  of  the  aspect  ratio 
effects. 

Several  factors  are  nonnegligible .  Manufacturing  techniques  are 
critical.  Surface  finish  on  the  model,  spanwise  twist,  and 
model-thickness  tolerances  produce  significant  effects  on  measured 
quantities  and  require  particular  attention  before  satisfactory  results 
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are  obtained.  The  boundary  layer  causes  a  horseshoe  vortex  to  fora  on 
each  airfoil  wall  intersection.  This  vortex  represents  a  departure  from 
two-dimensional  conditions  and  must  be  considered  (and  corrected)  before 
the  balance  can  give  an  accurate  force  representation. 

The  fourth  objective  concerned  the  way  to  implement  a  reliable 
force  balance  system.  The  recommended  way  is  to  modify  the  existing 
water  tunnel  balance  and  to  use  it  in  conjunction  with  Equation  (13) . 

Any  balance  which  measures  the  total  force  on  a  model  completely  span¬ 
ning  a  tunnel  with  an  incoming  boundary  layer  also  measures  an  addi¬ 
tional  increment  of  drag  due  to  the  three-dimensional  effect  of  the 
vortices  at  the  sidewall-hydrofoil  intersection.  The  use  of 
Equation  (13)  appears  to  correct  sufficiently  for  this  drag  increment. 
The  modifications  recommended  are  to  strengthen  the  weak  flexures,  to 
employ  only  two  force  cubes  as  sensors  with  the  cubes  being  oriented  so 
that  each  cube  measures  only  the  total  axial  or  total  normal  force,  and 
to  use  a  separate  torsion  element  to  measure  the  moment.  The  balance 
should  also  be  mounted  so  that  there  are  no  interactive  constants 
between  the  force  cubes. 

The  fifth  objective  was  to  formulate  a  data  reduction  process  to 
act  as  a  reference  procedure  for  future  testing  of  this  nature  at  ARL. 
The  data  reduction  process  is  outlined  in  Appendix  C.  Listed  in 
Appendix  D  are  tabulated  data  from  the  current  study.  An  error  analysis 
is  presented  in  Appendix  E. 

5.2  Recommendations  for  Future  Study 

Future  study  in  the  area  of  accurate  two-dimensional  force  measure¬ 
ment  should  include  further  testing  with  the  present  wind  tunnel 
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balance.  Several  different  airfoils  should  be  tested  to  establish  a 
large  data  base  for  establishing  the  correction  factor  (Equation  (jj)). 
Models  of  various  thicknesses,  various  aspect  ratios,  and  various 
degrees  of  camber  should  be  experimentally  investigated. 

Particular  attention  should  be  focused  on  the  thickness  variation 
in  follow-on  testing.  The  (t/c)^  term  of  Equation  (13)  requires  experi¬ 
mental  verification  because  the  three  studies  discussed  herein  used  only 
airfoils  of  12%  thickness.  Data  from  airfoils  of  differing  thicknesses 
were  not  available.  Since  Equation  (13)  is  empirical,  it  includes  the 
effects  of  drag  on  the  exposed  portion  of  the  force  transmitting  shaft 
(as  discussed  in  Section  3.4).  It  is  recommended  that  tare  drag  forces 
of  this  nature  be  further  investigated  to  determine  the  extent  of  their 
contribution  to  Equation  (13).  In  the  author's  opinion,  sufficient  data 
were  not  obtained  in  this  study  to  address  adequately  this  considera¬ 
tion. 

Another  recommendation  is  to  remove  the  horseshoe  vortex  physically 
by  a  suction  or  blowing  mechanism.  Data  taken  in  this  manner  should  at 
least  approach  the  two-dimensional  if  it  should  not  be  purely  two- 
dimensional  and  not  require  use  of  Equation  (13) .  A  comparison  test 
with  and  without  the  secondary  vortices  would  be  useful. 

It  is  also  recommended  that  the  pitching  moment  coefficient  be 
measured  by  some  procedure  other  than  by  differencing  the  measurements 
of  two  force  cubes.  The  addition  of  an  extra  force  cube  seems  to  exact 
a  penalty  in  drag  measurement.  Improving  the  flexure  setup  previously 
mentioned,  however,  could  possibly  overcome  this  difficulty.  Perhaps 
pitching  moment  measurements  could  be  made  with  a  torque  tube 
arrangement  or  even  with  a  different  balance  during  a  separate  set  of 
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experimental  runs.  After  these  recommended  tests  are  performed  in  the 
wind  ^unnel.  Equation  (13)  would  be  verified  (or  modified).  The  entire 
system  should  be  used  in  a  parallel  water  tunnel  test  program. 
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In  testing  of  this  nature,  the  presence  of  the  tunnel  walls  causes 
flow  conditions  in  the  tunnel  to  be  different  from  those  in  free  air. 
Conditions  are  altered  because  of  several  factors  and  the  usual  method 
of  correction  is  to  divide  the  interference  into  several  separate  com¬ 
ponents  and  simply  add  them.  Four  types  of  interference  are  considered 
here  as  traditional  corrections:  solid  blockage,  wake  blockage,  lift 
effect,  and  horizontal  buoyancy.  The  correction  terms  are  discussed  by 
Pope  (l4)  . 

Solid  blockage  refers  to  the  fact  that  the  presence  of  the  model  in 
the  test  section  reduces  the  area  through  which  the  fluid  must  flow. 
Because  of  the  Bernoulli  principle,  the  velocity  of  the  air  around  the 
model  is  increased  to  some  value  greater  than  that  which  the  model  would 
experience  in  free-stream  conditions.  Pope  refers  to  the  work  done  by 
Allen  and  Vincent!  (l5)  which  expresses  the  solid  blockage  term,  as 


£ 


sb 


A  o 


(Al) 


where 


wc.2 

0  “  (48)(h> 


A  = 


16 


f«l  "  p)U  +  £>) 1/2d  * 


x,y  *>  airfoil  coordinates;  c  *  airfoil  chord;  P  ■  no-camber  (basic) 
pressure  distribution;  and  h  »  tunnel  height.  Values  of  A  for  several 
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airfoil  families  are  plotted  as  a  function  of  thickness  in  Figure  6:8  of 

Reference  (l3) .  The  use  of  e  ^  will  be  discussed  later. 

Wake  blockage  is  also  a  continuity-related  condition.  In  the  wake 

behind  a  model,  the  velocity  is  lower  than  in  the  surrounding  stream. 

In  order  to  satisfy  conditions  of  continuity  in  a  closed  test  section, 

velocity  outside  the  wake  must  be  higher  than  the  velocity  far  upstream. 

In  a  free  fluid,  these  velocities  would  be  equal.  The  velocity  is 

therefore  accelerating  past  the  model  and  the  velocity  at  locations 

downstream  of  the  airfoil  is  higher  than  in  free-stream  conditions. 

Pope  gives  the  expression  for  wake  blockage,  e  ,  as 

w. 


(A2) 

/ 


where  t  »  (c/h)/4  and  C,  51  uncorrected  section  drag  coefficient.  The 

d 

u 

wake  blockage  and  solid  blockage  terms  represent  a  correction  to  axial 
velocity  past  the  airfoil  and  are  generally  combined  as 


e 


(A3) 


This  permits  a  velocity  correction  of  the  form 

V  -  Vu(l  +  e)  , 


(A4) 


where  the  subscript  u  represents  an  uncorrected  term  (velocity  in  this 
case).  The  Reynolds  number  then  becomes 

R  -  R  (1  +  e)  .  (A5) 

e  e 

u 
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Lift  effect  corrections  account  for  the  fact  that  the  tunnel  walls 
constrain  the  curved  streamlines.  The  effects  of  constraining  the 
streamlines  can  be  approximated  by  regarding  them  as  a  change  in  airfoil 
camber  and  angle  of  attack.  Pope  expresses  the  lift  effect  as  a  cor¬ 
rection  to  a,  C.,  C,,  and  C 

1  d  ml/4 

a  =  a  +  (--^--3-)-(C.  +  4C  )  (A6) 

u  2rr  *u  mlM 


C,  =  C*  (1  -  o  -  2e)  ,  (A7) 

u 


C 

ml/4 


C 

ml/4 


u 


oc* 

(1  -  2e)  +  , 


(A8) 


and 


Cd  •  Cd  11  -  3e,b  -  2S,b>  '  (A9> 

u 

where  C  *  pitching  moment  coefficient  about  the  1/4-chord  point. 
ml/4 

Horizontal  buoyancy  refers  to  the  additional  drag  created  by  the 
existence  of  a  static  pressure  gradient  in  the  tunnel  test  section.  As 
the  boundary  layer  builds  up  on  the  tunnel  walls,  the  axial  velocity 
increases  to  preserve  continuity  of  the  flow.  Along  with  this  increase 
in  axial  velocity  is  associated  a  decrease  in  static  pressure.  The 
downstream  pressure  is  lower  than  the  upstream  pressure  on  the  model; 
therefore,  a  spurious  drag  is  created  on  the  airfoil  which  must  be  sub¬ 
tracted  from  the  measured  drag.  Pope  gives  Allen  and  Vincenti's  expres¬ 
sion  for  this  buoyancy  drag,  D^,  as 
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where  dp/dfc  =  static  pressure  gradient  in  tunnel  without  an  airfoil. 
The  manner  in  which  these  corrections  are  incorporated  into  the  data 
reduction  process  is  discussed  in  Appendix  C. 


APPENDIX  B 
INTEGRATION  OF  WAKES 
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Drag  on  a  model  in  a  wind  tunnel  can  be  determined  by  comparing  the 
momentum  in  the  air  ahead  of  the  model  with  the  momentum  behind  the 
model.  The  drag  force  is  equal  to  the  loss  of  momentum  suffered  by  the 
air  as  it  passes  over  the  airfoil  or 

D  =  ^masS)  Change  in  velocity)  (Bl) 

sec 


or 


D  = 


pVda(V  -  V) 
o 


(B2) 


where  D  =  drag,  Vq  =  initial  airspeed  (ahead  of  airfoil),  V  *  final  air¬ 
speed  (behind  airfoil) ,  da  -  small  area  of  wake  perpendicular  to  air- 
stream,  and  p  ■  density.  Rearranging  terms  and  nondimensionalizing 
(details  in  Pope  (14)),  one  finds  that 


(  ( 


(B3) 


2  2 

where  q  ■  (l/2)pV  ,  qQ  *  (l/2)pVQ  ,  y  =  unit  area  term  measured  perpen¬ 
dicular  to  the  airfoil,  and  c  »  chord.  In  Equation  (B3) ,  the  term  in 
parenthesis  can  be  rewritten  as 


(B4) 
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The  drag  coefficient  can  now  be  determined  by  measuring  a  reference 
velocity,  Vq,  and  a  series  of  local  velocities  which  are  measured  by 
traversing  a  probe  through  the  airfoil  wake. 

The  wake  measurement  tests  were  set  up  as  shown  in  Figure  Bl.  The 
traversing  probe  was  a  kiel  probe  which  accurately  measures  dynamic 
pressures  for  relative  wind  angles  within  ±35°  of  the  probe  centerline. 
The  static  pressure  reading  taken  at  the  reference  probe  location  was 
used  for  both  V  and  Vq  computations.  This  is  acceptable  since  both 
probes  are  located  the  same  distance  downstream  of  the  panel  leading 
edges.  It  proved  to  be  very  difficult  to  make  local  static  pressure 
readings  from  the  traversing  probe  due  to  streamline  curvature  in  the 
airfoil  wake. 

The  movable  probe  measured  local  velocities  (dynamic  pressures)  as 
it  was  traversed  in  0.1-in.  increments  vertically  through  the  wake. 

This  profile  of  velocity  ratios  (such  as  depicted  in  Figure  14)  was  then 
integrated  using  the  trapezoidal  rule  to  give  the  value  of  the  double 
integral  in  Equation  (B3) .  The  computer  program  and  the  associated 
variables  are  listed  in  Appendix  C. 
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Figure  Bl.  Wake  traverse  installation. 


APPENDIX  C 


COMPUTER  PROGRAMS  AND  DATA  REDUCTION  PROCEDURE 
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The  two  computer  programs  which  were  used  to  calculate  the  force 
coefficients  and  integrate  the  airfoil  wakes  are  presented  herein.  The 
force  coefficient  program  takes  measured  temperatures,  pressures,  angle 
of  attack,  and  balance  voltages  and  determines  corrected  C^,  C^,  and  Re. 

The  final  wind  tunnel  balance  configuration  is  represented  in  the 
top  of  Figure  Cl.  Tension  members  A  and  B  are  calibrated  (KA  ■  cali¬ 
bration  constant)  such  that  the  tensile  force  and  direction  along  the 
ribbon's  long  axis  can  be  determined,  e.g., 

F°rCeA,B(lb)  "  KA,B(^U)V°lt8A,B  *  (C1) 

where  the  sign  of  A's  or  B's  voltage  reading  would  indicate  whether  the 
force  was  directed  up  or  down  along  the  ribbon  (indicated  by  double¬ 
headed  arrows  in  Figure  Cl.  In  order  to  determine  lift  and  drag,  nor¬ 
mal  and  axial  forces  must  be  transformed.  The  lower  portion  of 
Figure  Cl  shows  the  appropriate  vector  diagram.  To  produce  a  lift  force 
or  vector,  £,  ribbons  A  and  B  must  be  stressed  as  vectors  and  as 
shown.  Likewise,  to  produce  a  drag  vector,  5,  ribbons  A  and  B  must  be 

stressed  as  vectors  D,  and  5_.  Therefore,  when  both  lift  and  drag  exist 

A  B 

on  an  airfoil,  ribbon  A  will  reflect  the  vector  sum  of  and  and 
ribbon  B  will  reflect  the  vector  sum  of  and  5_.  To  determine  how 

O  D 

much  of  each  ribbon's  total  measured  force  contributes  to  lift  or  drag, 
two  simultaneous  equations  must  be  solved.  From  the  vector  diagram  of 
Figure  Cl,  it  can  be  shown  that 

D(lb)  =  F.(lb)sin  o  +  F  (lb)cos  a  (C2) 

A  D 


and 
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% 


Figure  Cl.  Schematic  diagram  of  relationship  between  aerodynamic 
forces  and  cantilever  balance  elements.  Ribbon  B  is 
aligned  with  the  airfoil  chordline  and  rotates  with 
the  airfoil.  Ribbon  A  is  perpendicular  to  ribbon  B. 
In  the  lower  portion  of  the  figure,  the  circled 
values  represent  sign  conventions  for  the  individual 
ribbons. 
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L (lb)  =■  FA(lb)cos  a  -  FB(lb)sin  a  ,  (C3) 


where  F  represents  the  forces  measured  by  the  two  ribbons  (determined 

Aj  O 

as  per  Equation  (Cl) . 

The  force  coefficient  program  then  follows  this  sequence: 

a.  Read  barometric  pressure,  number  of  events,  transducer  con¬ 
stant,  airfoil  chord,  and  static  pressure  gradient. 

b.  Calculate  repeatable  terms:  o,  £  .  ,  t. 

SD 

c.  Establish  elements  of  simultaneous  equations  subroutine. 

The  matrix  is  of  the  form 


»1 

xi 

lBJ 

A2,1A2,2 

N 

(C4) 


where  A.  ■  trignometric  functions  of  a  in  Equations  (C2) 

1  tJ 

and  (C3);  ■  forces  A  and  B  in  pounds  as  determined  by 

Equation  (Cl);  and  -  unknown  forces,  lift  and  drag. 

d.  Determine  uncorrected  lift  and  drag  in  pounds. 

e.  Determine  buoyancy  drag  increment  from  Equation  (A10), 


_  6h  ,  dp 

db - Ao 


Ao  .. 
n  di 


f.  Nondimens ionalize  coefficients 


C*  " 


u  (l/2)pV  sc 


(C5) 


(C6) 


and 
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D-D 


B 


(1/2) pV  sc 


(C  7 ) 


g.  Correct  coefficients  and  a  for  solid  blockage,  wake  block¬ 
age,  and  lift  effect: 


a 


a  + 

u  2tt  il 


u 


^Cm  1/A  ) 
u 


(C8) 


(The  moment  term  of  Equation  (C8)  was  neglected  in  this 
study  since  moment  forces  were  not  measured.  It  is  in¬ 
cluded  here  in  the  event  that  moments  are  measured  in 
future  testing.) 


C,  -  C,  <1  -  »  -  2c) 

u 


(C9) 


and 


C.  (1  -  3e  .  -  2e  ,  ) 
d  sb  wb 

u 


(CIO) 


Re  "  Re  (1  +  esb  +  ewb)  ‘  (C11> 

u 

The  program  is  listed  with  an  explanation  of  variables  at  the  end 
of  the  appendix. 

The  wake  integration  program  is  listed  here  alone.  The  wake  inte¬ 
gration  program  is  discussed  in  Appendix  B. 
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VARIABLES  IN  FORCE  COEFFICIENT  PROGRAM 

A  =  matrix  coefficient  (Equation  (C4)) 

ALFA  =  airfoil  angle  of  attack  (in  degrees) 

ALFAC  =  airfoil  angle  of  attack  with  classic  corrections  applied  (in 
degrees) 

ANGL  =  angle  formed  by  top  edge  of  balance  and  horizontal;  equal  to 
ALFA  in  final  balance  configuration  (in  degrees) 

B  =  matrix  coefficient  (Equation  (C4 ) ) 

BAR  =  barometric  pressure  (in  inches  of  mercury) 

BETA  =  ALFA  (in  radians) 

C  =  airfoil  chord  (in  inches) 

CD  =  drag  coefficient 

CDC  =  drag  coefficient  with  classic  corrections  applied 

CL  =  lift  coefficient 

CLC  =  lift  coefficient  with  classic  corrections  applied 

CON  =  transducer  constant  (in  inches  of  water/volt) 

D  =  drag  (in  pounds) 

DB  =  buoyancy  drag  (in  pounds) 

2 

DPDL  =  static  pressure  gradient  ( - 

(l/2)pV^ 

ESB  ■  solid  blockage 

EVENT  *  test  number 

L  ■  lift  (in  pounds) 

N  ■  number  of  angle-of-attack  positions  surveyed 

PI  -  3.14159 


R 


gas  constant 
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RDELP  =  reference  velocity  measurement,  total  pressure-static  pressure 

RRE  =  reference  Reynolds  number 

REC  =  reference  Reynolds  number  corrected 

RRHO  =  reference  density 

RVEL  =  reference  velocity 

RXNU  =  reference  dynamic  viscosity 

SIGMA  =  o  as  in  Equation  (Al) 

TAU  =  t  as  in  Equation  (A2) 

TEMP  =  temperature  (in  degrees  Fahrenheit) 

VA  =  ribbon  A's  voltage  (in  volts) 

VB  *  ribbon  B's  voltage  (in  volts) 

XKA  *  ribbon  A's  calibration  constant  (in  pounds/volt) 

XKB  *  ribbon  B's  calibration  constant  (in  pounds/volt) 


FORCE  COEFFICIENT  PROGRAM 


PROGRAM  FORC £  2 
REAL  1(501 
INTEGER  EVENT 

01  ME  NS  I  ON  TEMPI  SO) .ALF A( 50) *VA( 50  I • V3<  SOI «0( 50) *ANGL1 50) * 
l ALFACI 50) .CL ( 501 «C0(50I .OBI  501 .CLC  <  501 «COC) 50  >  *REC I  50) .RRMOI SO) . 

2R0ELPI 50) »RVEL(  50) .RXNUI SOI .R«E ( 50) .BET  AC  50  I .  Al 50) .B( 50 l 
RE  AO (A, 99)  N.EVENT.SaR  .CON.C  .DPDL 
99  F0RMAT(2I5.4F10.5> 

PI  =  3.1415927 

SIGMA=( P1992/48. >9(C/48. 1**2 
ESB  -  ( . 24 ! * S 1 GMA 
TAU  =  (C/43.1/4. 

R  *  1715. 

00  100  I =1 ,N 

REA0(4, 101  I  ALFA<  I  1  «  ANGL  (  ll.VAUl.VBdl.RDELPUl.TFMPUl 
LOl  FORmaT ( 6F 1 0. 5 ) 

TEmp(I)  =  (TEMPI  11*100. 0)»460.0 
100  CONTINUE 

WR I T  E (  3.5981 

593  FORMAT ( 'G'.SX.'ALFA'.SX.'ALFAC'.bX.'CL'.ax.'CLC'.TX.'CD'.RX.'CDC'. 
17X,'REC .SX.'ULBI • ,5X,'D(LB>' , 3X ,  •  RVEU  F/S  I • ) 

XKA  =  -.17598 
XK8  =  -.14485 
VAVE  *  0.0 
00  10  1=1. N 

BETA ( I  I  =  (ANGL! I ) /  1 7 . 777 7 7 B ) * ( P 1/ 1 BO.O  ) 

All)  r  SINIBETAI  I  I ) 

A( 2 )  =  COSIBETAI 1 ) I 
A (  31  =  COS(BETA(  I  1  I 
A( 4)  =  SINIBETAI I ))*) -1.0) 

BUI  =  V  A  (  I  I  /XKA 
8(21  =  VBIII/XKB 
CAU  S1M3I  A.B.2.KS) 

IFIKS.EQ.I )  MRITEI5.597) 

597  FORMAT) 'MATRIX  ERROR') 

DID  *  B(l) 

U  II  =  B  (  2  I 

ROELPI I  I  =  ROELP( f )9C0N«5. 20128 

RRMOI  1 1  =  .O02378*IBAR/29.92  )<■)  519.0/TEMPI  1  )  I 

RVEL(I)  =  SORT! (2.SRDELPI I  I l/RRHOI I) I 

RXNUII)  =  ( • 3Tt-6 ) /RRMOI I  I 

RRE(I)  =  (RVELII )«IC/12. ) l/RXNUl I) 

Ofil  I  1  =  lb.*(48./12.  ) «*2. /PI  )9.24*SIGMA*(0P0l*.5»RRHO(  I  )*RV6L«  I 
12.4(0/12. )  I 
VAVE  *  VAVE'RVEMII 
10  CONTINUE 

VAVE  *  VAVE/N 
00  11  I=i.N 

CUII  *  U  1 )/(  .59RRH0I  1 1 « (  VAVE**2.  I  «1 . 53646*1  C/12.  I) 

COU)  =  tOm-t)B(II  I/(.5*RRH0(I)*VAVE**2.*I.53646*(C/I2.)I 

REC(I)  *  RRE 1  I ) *1 1 . ♦ ( ESB»TAU*C0< I  I  I  I 

CICIII  *  CUII*U. -SIGMA-12. *ES8)-I2.*TAU*C0|I>>) 


CtKUl  =  CDI  1 1*1  l.-<3.*ESBI-(2.eTAUeC0<II  II 
Al FA ( 1  I  s  ALFA! I  1/17.777776 

ALFACUI  =  ALFA  (  l  )  ♦  1  I57.29CSI  GMA/2 ./P I  I  »CL  (  I  )  I 

WRITE! 3, 105)  ALFA!  1 1  ,ALFaC<  1 1  .CU  1 1  ,CLC  I  11  .CO  ( I)  .CDC  l  II  .REC  I  I )  • 
lLIII.OIII,RVtL<ll 
105  FORMAT! *  0  * .6F 10 .5 , F 10. I « 3F 10. 5  I 
11  CONTINUE 

WR1TEI3.599I  EVENT, VAVE 

599  FORMAT ( *  O' , • RUN  NO. • , 2X . 1 5 , 5* , • AVC  VEL  I S • , 2X ,F 10. 5  I 
STOP 
END 
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VARIABLES  IN  WAKE  INTEGRATION  PROGRAM 

AREA  =  area  under  velocity  profile  curve  being  integrated 

BAR  =  barometric  pressure  (in  inches  of  mercury) 

C  =  airfoil  chord  (in  inches) 

CD  a  two-dimensional  drag  coefficient 

CON  =  transducer  constant  (in  inches  of  water/volt) 

DELP  =  local  velocity  measurement;  local  total  pressure-reference 
static  pressure  (in  volts) 

DLPF  =  reference  velocity  measurement;  reference  total  pressure- 
reference  static  pressure  (in  volts) 

END  =  integration  limit 
EVENT  =  test  number 

IS  ■  sign  of  voltage  measurement  (+  or  -) 

ISZER  =  amplifier  offset  voltage 

N  =  number  of  positions  in  wake  at  which  measurements  were  taken 

P  =  absolute  static  pressure  (in  pounds/cubic  foot) 

PDT  *  local  total  pressure  (in  volts) 

PIS  =  reference  static  pressure  (in  volts) 

« 

PIT  =  reference  total  pressure  (in  volts) 

R  *  gas  constant 

RATIO  *  ratio  of  local  velocity  to  reference  velocity 

RHO  *  density 

START  °  integration  limit 

STAT  ■  static  pressure  difference  from  atmospheric  (in  pounds/cubic 
foot) 

■  temperature  (in  degrees  Fahrenheit) 


TEMP 
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VEL  *  local  velocity 
VELIF  =  reference  velocity 

VINTGL  *  rectangular  Increment  of  wake  velocity  curve  which  is  summed  to 
measure  area  under  wake  velocity  curve 
VOLT  **  any  of  the  voltage  measurements  as  recorded  on  paper  tape. 

This  term  is  used  for  subtracting  out  instrument  offsets  from 
all  measured  quantities. 

ZERVT  =  pressure  transducer  offset  voltage 
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WAKE  INTEGRATION  PROGRAM 


PROGRAM  WARE 
INTEGER  START,  EVE..T 

01  MENS  I  UN  PIT!  1501  .PISUfo  I  .POT!  150.  TFmP(  ISO) 

DIMENSION  RAT IP! ISO) ,STAT I  ISO) »  P  I  ISO  I ,0£LPI l  5G I »  PhD  ! ISO! 

OIMENSION  V INTGL I ISOI , VEL l 1 50  I . VOL  T I  2001 . IS! lol 
DIMENSION  OLPr I 20J I .VEL 1= I 2uO I 
RE  ADI  4. SO  IN , F VENT, BAP, CON, C  .START,  I  END 
50  FORMAT I2I5,3F10.5,2IS| 

WRITE(3.88«I 

888  FORMAT!///, ■  • , 5* , • P SN • , 7x , • R A T 1 0 • . 1 C A . • VE L • , 1 0* , • VREF • I 
DO  100  1=1. N 

RE  ADI  101  I  IS2ER.2ERVT.  I  I  SI  J)  .VOL  T ( Jl ,J* l ,81 

10 1  FORMAT! 8X« 5 ( I l,F8.3,5X I .1 I.F8.3I 
CALL  CONVRT  (  8  ,  I  S  .  VOL  T  ,  I  lP  V  T  .  I S7E»  I 
00  102  J =2.4 

VJLTI JI=VOLT( JI-VOLTI 1  I 

102  CONTINUE 

PIT! II  =  VOL Tl 2  I 
POT! I )  =  VDLT! 3  I 
PISI I  I  =  VOL T ( 4  I 

TEMP! I ) *1  VOLT! 5  IS  100.0  I *460.0 
100  CONTINUE 
R  =  1 7 1 5.0 

C  BAROMETRIC  PR  NOW  CONVTO  FROM  IN.  m£RC  10  LB/SO  FT 

3AR=3ARS70.7328 
DO  3  I=l»N 

C  STATIC  PR  NOW  CQNVTD  FROM  VOLTS  TO  LP/SW  FT 

STAT! 1 1 =P I S ( I tsC0N«s.20128 
PI  I  I *BAR»  ST  AT  I  I  I 

C  ABSOLUTE  STATIC  PR  IS  NOW  IN  LB/SO  FT 

DEEP  I  11  =  1  POT  I  I )-° I  SI  I  I ) *C0N  =  5 • 20 1 26 
DLPF( 1 1 =( PIT! 1 1 -PI  SI  I  I l*CONS5.20128 
RHOt I  1  =  1  PI  I l/IRSTEMPI I >1)0  32. 174 
VEL I  I  I  =  SORT ( I 2.00EEP! I  I o 32 . 1 74 ) /RhO I  I  I  I 
VELIFI 11  =  SORT  I ( 2.S0EPF 1 1 |S32< 1 74  1/RMOI 1 1  I 
RAT  10 1 1  I =VEL I l I /VELIFI I ) 

WRITEI3.510) I .RATIOI I  I «VFEI 1 • «VEL  I  e  I I  I 
510  FORMA  T I •  *»5X.I3.F12.4.2X.F12.A*2X,F12.4> 

3  CONTINUE 

C  CONVERT  BAR  PR  8ACX  INTO  INCHES  FOR  PRINTOUT 

BAR=9AR/70.7325 
AREA=0.0 

00  70  I=START.I ENO 

VINTGLI I) »IH RATIO! II -RATIOI I  I »S2 I ♦ I R AT  1 01  I  *  1 1 -RAT  10! I ♦ 1 1 **2 >1 / 

12. IP! .1/12.1 
AREA=AREA*VINTGU 1 1 
70  CONTINUE 

CD* 1 2./I C/12. ) I  PARE A 
WRITE  I  3.501  ICO 

501  FORMAT  I *0* » 8X« *0P AG  COEFF  I  S’ . 2X.F8.4/I 
WRITE  1 3 1 108  I  BAR .EVENT 

108  FORMAT!*  *  *8X , *8AR0M  PRESS  I S * . 2X , F6. 3 , 1 X , » IN.  MG *. 5X, • SURVEY  NO*. 
12X.I5/I 
STOP 
ENO 
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DATA  TABLES 
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The  tabulated  data  included  here  consist  of  lift  and  drag  data 
measured  by  the  doubly  supported  balance,  lift  and  drag  measured  by  the 
cantilever  balance  with  both  the  9-  and  11-inch  chord  airfoils  and  with 
and  without  disk,  and  mid-span  wake  traverses  which  produced  values 
for  various  a's. 
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TABLE  D1 

Lift  and  Drag  Data  -  Doubly  Supported  Balance.* 


V  =  110  f/s 
Re  =  330,000 
c  =  6  in. 


a  (deg) 

CL  (±0.0024) 

CD  (±0.0002) 

8.26 

0.773 

0.0170 

6.47 

0.626 

0.0136 

4.91 

0.480 

0.0114 

2.70 

0.215 

0.0087 

1.04 

0.053 

0.0084 

*data  corrected  for  tunnel  interference  effects  and 
Equations  (1)  and  (5) 
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TABLE  D2 

Lift  and  Drag  Data  -  Cantilever  Balance  - 
9-in. -Chord  Airfoil,  11-in.  Disk.* 


V  =  70  f/s 
Re  =  330,000 
c  =  9  in. 


a  (deg) 

CL  (±0.0020) 

CD  (±0.0006) 

8.77 

0.794 

0.0271 

6.92 

0.682 

0.0219 

5.18 

0.549 

0.0174 

2.90 

0.338 

0.0131 

0.91 

0.086 

0.0101 

-1.11 

-0.086 

0.0102 

-3.17 

-0.346 

0.0122 

-5.19 

-0.530 

0.0151 

-7.16 

-0.686 

0.0194 

-9.37 

-0.805 

0.0265 

*data  corrected  except  for  disk  tare  forces  and 
Equation  (13) 
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TABLE  D3 

Lift  and  Drag  Data  -  Cantilever  Balance 
(Disk  Tare  Forces) .* 


V  =70  i/s 

Re  =  330,000 

disk  dia  =  11  in. 


a  (deg) 

CL  (±0.0014) 

CQ  (±0.0004) 

9.34 

0.0094 

0.0021 

6.26 

0.0078 

0.0014 

3.21 

0.0055 

0.0018 

0.47 

0.0039 

0.0012 

-3.43 

0.0027 

0.0009 

-6.17 

-0.0008 

0.0014 

-9.01 

-0.0004 

0.0020 

Coefficients  nondimensionalized  by  c  =  9  in. 
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TABLE  LA 


Lift  and  Drag  Data  -  Cantilever  Balance  - 
18-in. -Chord  Airfoil,  No  Disk.* 


V  =  37  f/s 
Re  =  330,000 


a  (deg) 

C  (±0.0035) 

CD  (±0.0007) 

-10.38 

-0.838 

0.0289 

-8.55 

-0.736 

0.0207 

-6.30 

-0.593 

0.0150 

-A.  12 

-O.A27 

0.0113 

-1.98 

-0.17/ 

0.0087 

-0.26 

-0.018 

0.0076 

2.08 

0.195 

0.0096 

3.10 

0.32A 

0.0102 

3.95 

0.A15 

0.0121 

6.03 

0.575 

0.0158 

6.17 

0.595 

0.0158 

8.02 

0.715 

0.0220 

9.44 

0.799 

0.0276 

10.29 

0.830 

0.0329 

*all  corrections  applied  except  Equation  (13) 
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TABLE  D5 

Lift  and  Drag  Data  -  Cantilever  Balance  - 
9-in. -Chord  Airfoil,  No  Disk.* 


V  =  70  f/s 
Re  =  330,000 


a  (deg) 

CL  (±0.0038) 

CD  (±0.0005) 

11.41 

0.909 

0.0460 

10.78 

0.890 

0.0400 

10.44 

0.875 

0.0371 

8.78 

0.790 

0.0273 

8.12 

0.757 

0.0249 

6.75 

0.672 

0 . 0206 

6.01 

0.613 

0.0179 

5.17 

0.559 

0.0163 

3.95 

0.463 

0.0138 

2.98 

0.358 

0.0120 

-0.51 

-0.026 

0.0089 

-3.34 

-0.361 

0.0109 

-6.34 

-0.618 

0.0150 

-9.09 

-0.796 

0.0233 

*all  corrections  applied  except  Equation  (13) 
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WAKE  TRAVERSE  DATA 

The  following  computer  printouts  are  wake  traverse  data  for  the 
9-in. -chord  airfoil  at  various  angles  of  attack.  "Position"  indicates 
vertical  location;  "VEL"  is  local  velocity.  "VREF"  is  the  reference 
velocity  measured  at  the  same  instant  as  the  local  velocity  to  eliminate 
any  temperature  change  effects  during  the  test.  "RATIO"  is  VEL/VREF  and 
is  the  value  (V/V^)  used  in  Equation  (B2).  Listed  in  Table  D6  are  the 
starting  positions  of  each  traverse  (Position  1) .  Positions  2-51  move 
toward  the  tunnel  ceiling  in  0.1-in.  increments.  All  velocities  are  in 
feet/second.  Tunnel  height  is  48  in.  Re  =  330,000.  The  probes  are 
located  3.1  c  behind  the  airfoil. 
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TABLE  D6a 


Survey  45,  a  =  0.3°. 


PSN 

RA1  10 

VFL 

V«£F 

1 

1 .0023 

67. 8981 

67.6930 

2 

1.0034 

67.8666 

67.6338 

3 

0.9977 

67.9Q59 

6e .0600 

4 

1.0017 

68.0405 

67.9244 

5 

1 .0006 

68.0018 

67.9631 

6 

0.9994 

68.0533 

68.0920 

7 

1 .0000 

67.9495 

67.9495 

d 

0.9994 

68.0721 

68.1108 

9 

1.0011 

68.1620 

68.0847 

10 

1.0017 

68.2780 

68.1622 

1 1 

1 .0023 

68.3297 

68 .1753 

12 

1 .0006 

68.0585 

68.0198 

13 

1 .0017 

68.1361 

68.0200 

14 

1.0000 

68.0586 

68.0586* 

15 

0.9966 

67.5655 

67.7992 

16 

0.9954 

67.4544 

67.7664 

17 

0.9931 

67.2318 

67.700B 

Id 

0.9868 

67.1934 

68.0901 

19 

0.9814 

66.4883 

67.7459 

20 

0.9767 

66.1243 

67.7008 

21 

0.9715 

65.9032 

67 • R  367 

22 

0.9661 

65.4953 

67.7914 

23 

0.9639 

65.5018 

67.9535 

24 

0.9587 

65.7164 

68.5676 

25 

0.9662 

65.5564 

67.8517 

26 

0.9697 

65.8034 

67.9577 

27 

0.9744 

66.0833 

67.8185 

28 

0.9779 

66. 3286 

67.8246 

29 

0.9884 

66.9173 

67.7018 

30 

0.9919 

67.2778 

67.8251 

31 

0.9936 

67.2836 

67.7140 

32 

0.9983 

67.6299 

67.7469 

33 

0.9988 

67.5124 

67.5905 

34 

1.0012 

67.6028 

67.5246  * 

35 

1.0012 

67.6878 

67.6097 

36 

1.0035 

67.7657 

67.5313 

37 

1.0012 

67.6093 

67.5311 

38 

1 .0000 

67.7266 

67.7266 

39 

1.0012 

67.6939 

67.6158 

40 

1.0023 

67.7383 

67.5826 

41 

1.0023 

67.8623 

67.7062 

42 

1.9006 

67.4712 

67.4320 

43 

1.0012 

67.6279 

67.5497 

44 

1.0041 

67.7061 

67.4323 

45 

1.0029 

67.7517 

67.5563 

46 

0.9971 

67.4838 

67.6794 

47 

1.0000 

67.5227 

67.5227 

48 

1.0023 

67.6733 

67.5168 

49 

1.0006 

67.7191 

67.6800 

SO 

1.0000 

67.5621 

67.5621 

51 

0.9994 

67.7313 

67.7703 

DRAG 

COEFf  IS 

0.0075 

BARQM 

PRESS  IS 

28.960  IN.  HG 

SURVEY  NO 
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TABLE  D6b 


Survey  46,  a  =  1.57°. 


PSN 

RATIO 

VEl 

V«£F 

1 

1  .0196 

67.3965 

66.0989 

2 

1.0023 

67 •*704 

67.4150 

3 

1 .0035 

67 .7C57 

67.4727 

4 

1 .00  l  1 

67.8350 

67.7575 

5 

1.0023 

07.9673 

67.8026 

b 

0.9960 

67.6601 

67.9315 

7 

l .0011 

67.7886 

67.7112  * 

8 

0 .990  3 

67.8010 

67.9173 

9 

0  .9983 

67.8077 

67.9240 

10 

0.9994 

0 7 • 8654 

67.9042 

11 

0.9960 

©7 ,5091 

67. 7yl4 

12 

0.9948 

67  .5219 

67.8719 

13 

0.992b 

67.4108 

b7  .9167 

14 

0.9351 

66.9590 

67.9742 

15 

0.9838 

66.5975 

67.6958 

16 

0.9754 

66 .5706 

68.2513 

IT 

0 ,9o9  3 

65.9470 

68.0323 

18 

0.9b69 

65.7121 

b7 .9603 

19 

0.9623 

65.4774 

68 .0446 

20 

0.9571 

65.2289 

68.1540 

21 

0.9545 

64.9101 

68.0047 

22 

0.95e3 

65.3626 

68.2058 

23 

0.9571 

65.2883 

68.2125 

24 

0.9647 

08. 6901 

bB .0956 

25 

0.9741 

66.3419 

6e.lOB3 

26 

0.9787 

66. 5404 

67.9918 

27 

0.9880 

67.2631 

68.0822 

23 

0.9925 

67.2233 

67.7318 

29 

0.9943 

67.8488 

68.2367 

30 

1 .0000 

68.1332 

68.1332  * 

31 

0.9994 

67.7437 

67.7827 

32 

1.0006 

67.7825 

67.7436 

33 

1.0023 

67.9901 

67.8345 

34 

0.9994 

67.8738 

67.9127 

35 

1.0017 

68.0228 

67.9063 

36 

1.0035 

67.8739 

67.6400 

37 

1.0017 

67.8470 

67.7300 

38 

1.0023 

68.0808 

67.9254 

39 

1 .0034 

67.9702 

67.7365 

40 

1.0029 

68.0804 

67.8860 

41 

1.0040 

68.2417 

67.9701 

42 

1 .0006 

67.9381 

67.8992 

43 

1.0023 

68.0928 

67.9373 

44 

1.0011 

60.0607 

67.9830 

45 

l .0006 

67.8660 

67.8270 

46 

1.0017 

68.1833 

68.0669 

47 

0.9994 

67.9885 

69.0274 

46 

0.9989 

68.1441 

68.2217 

49 

I .0046 

68.1441 

67.8329 

50 

1.0000 

67.7998 

67.7998 

51 

1.0017 

68.4291 

68.3130 

DRAG 

C06FF  IS 

0.0099 

BAR  0“ 

PRESS  IS 

28.960  IN.  HG 

SURVEY  N 

/ 


46 
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TABLE  D6c 

Survey  47,  a  =  4.00°. 


PSN 

RATIO 

VEl 

VREP 

1 

1.0274 

67.2530 

65.7140 

2 

1 .0029 

67.6356 

67.4411 

3 

1 .0006 

67.5309 

b  7 .4920 

4 

1.0006 

67.6993 

67. 6605 

5 

1 .0000 

67.7573 

67.7573 

6 

1.0023 

67.8017 

67.6464 

7 

1.0011 

67. 8984 

67.8208  * 

8 

0.9954 

67.5550 

67.8660 

9 

0.9971 

67.4892 

67.6840 

10 

0.9819 

66.2346 

67.4558 

U 

0.9743 

65.6419 

67.3060 

12 

0.9640 

65.0026 

67.4295 

13 

0.9508 

64.1640 

67.4874 

14 

0.944 1 

63.4199 

67.1738 

15 

0.9425 

63.4317 

67.3040 

lb 

0.9424 

63.3904 

67.2651 

17 

0.9454 

63.8936 

67.5843 

18 

0.9604 

64.8418 

67.5123 

19 

0.9729 

65.5761 

67.4014 

20 

0  . 9d03 

66.3027 

67.6032 

21 

0.9918 

66.7778 

67.3287 

22 

0.9977 

67.1382 

67.2983 

23 

1.0006 

67.4649 

67 .42  58* 

24 

1.0006 

67.6994 

67.6604 

25 

1.0035 

67.6992 

67.4648 

26 

1 .0017 

67.7510 

67.6340 

27 

1.0029 

67.6728 

67.4773 

28 

1.0029 

67. 7964 

67.6013 

29 

1.0012 

67.6072 

67.5290 

30 

1.0041 

67.7758 

67.5022 

31 

0.9983 

67.4899 

67.6072 

32 

1.0029 

67.7698 

67.5745 

33 

1.0017 

67.6851 

67.5680 

34 

1.0012 

67.6204 

67.5421 

35 

0.9983 

67.5092 

67.6266 

36 

0.9994 

67.3581 

67.3974 

37 

1.0029 

67.6263 

67.4305 

38 

1  .0000 

67.6387 

67.6387 

39 

1.0052 

67.8729 

67.5213 

40 

1.0035 

67.7557 

67.5211 

41 

1.0035 

67.7620 

67.5274 

42 

1.0017 

67.7164 

67.5992 

43 

1.0012 

67.4876 

67.4092 

44 

1.0000 

67.7623 

67.7623 

45 

1.0023 

67.7417 

67.5852 

46 

1.0000 

67.6575 

67.6575 

47 

1.0023 

67.8136 

67.6574 

48 

1.0052 

67.9242 

67.5728 

49 

1.0023 

67.6695 

67.5128 

50 

1.0017 

67.7086 

67.5912 

51 

1.0000 

26.6953 

26.6953 

DRAG 

COEFF  IS 

0.0097 

baron 

PRESS  IS 

28.940  IN.  HG 

SURVEY  NO 
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TABLE  D6d 

Survey  48,  a  =  5.92°. 


PSN 

RATIO 

vet 

1 

1.0114 

68.8497 

2 

0. 9988 

67.2359 

3 

1 .0029 

67.4379 

4 

l .0000 

97.3263 

5 

0.9988 

67.2996 

6 

0.9925 

67.1098 

7 

0.9878 

66.7213 

8 

0.9749 

65.9762 

9 

0.9677 

65.4994 

10 

0.9551 

64.6051 

1  1 

0.9378 

63.3761 

12 

0.9287 

62.8306 

13 

0.9322 

63.0049 

141 

0.9353 

63. 1727 

15 

0.9423 

63.8157 

1 6 

0.9551 

64.6694 

17 

0.9621 

65.4486 

18 

0.9/45 

66.2990 

19 

0.9872 

66.8219 

20 

0.9954 

67. 4469 

21 

1.0000 

67.8053 

22 

1.0029 

68.1300 

23 

0.9988 

67.7784 

24 

1.0029 

68.0970 

25 

l .0006 

67.9740 

26 

0.9966 

67.7852 

27 

1.0023 

6". 0256 

28 

1 .0023 

67.7912 

29 

1.0006 

67.7518 

30 

1.0029 

67.7975 

31 

1.0017 

68.C379 

32 

0.9977 

67.4108 

33 

1.0040 

67.9142 

34 

1.0029 

67.7309 

35 

1.0017 

67.6919 

36 

1.0012 

67.6197 

37 

1.0006 

67.7828 

38 

1.0058 

67.7825 

39 

1.0035 

67.9779 

40 

1 .0000 

67.7825 

41 

1.0029 

67.9778 

42 

1.0029 

67.9391 

43 

1.0006 

67.7500 

44 

1.0000 

67.8343 

45 

0.9977 

67.6113 

46 

1.0006 

67.7498 

47 

1.0017 

67.9584 

46 

1.0017 

67.6776 

49 

1.0000 

67.6B37 

50 

0.9965 

67.4870 

51 

1.0029 

67.8013 

ORAG 

COEFF  IS 

0.0120 

BARON 

PRESS  IS 

28.940  (N. 

"WAT  /■  — '  -fr 


VREF 

66.0957 

67.3144 

67.2418 

67.3263*  1 

67.3791 

67.6198  1 

67.5483 

67.6778 

67.6840  * 

67.6451 

67.6787 

67. 6566  I 

67.5848 

67.5456 

67.7210 

67 .7086 

69.0271 

68.0336 

67.6882 

67.7603 

67.8053* 

67.9352 

67.8565 

67.9022 

67.9350 

68.0193 

67.8696 

67.6347 

67.7127 

67.6017 

67.9209 

67.5679 

67.6404 

67.5349 

67.5743 

67.5412 

67.7436 

67.3900 

67.7433 

67.7825 

67.7824 

67.7436 

67.7108 

67.8343 

67.7681 

67.7106 

67.8412 

67.5598 

67.6837 

67.7226 

67.6053 


survey  no 


48 


128 


TABLE  D6e 


Survey  49,  a  *  9.95°. 


PSN 

RATIO 

VEC 

VR£f 

1 

1 .0018 

66.9651 

66,8670 

2 

uooao 

67.0167 

67.0167 

3 

1 .Q0Z9 

67.1530 

66  «  956*. 

4 

1  .0013 

66 • 9o8 1 

66.9698 

5 

1.0000 

67.25C2 

67.2502 

6 

1.0035 

67.3802 

67.1667 

r 

1 .0012 

67.2669 

67.1903 

e 

1.0061 

67.31 38 

67.0386 

9 

1.0012 

67.3202 

67.2617 

to 

0.9988 

67.1359 

67.2165 

u 

1.0000 

67.3780 

67,3780  * 

12 

0.9988 

67.6261 

67.5026 

13 

0.9912 

66.9000 

67.3910 

u 

0.9866 

66.5677 

67.6751 

15 

0.9813 

66. 26 1 6 

67,5269 

16 

0.9765 

65.8615 

67.6697 

17 

0.9657 

65.1612 

67.6961 

18 

0.9565 

66.6172 

b7 • 36°5 

19 

0.9661 

63.8607 

67.6629 

20 

0.9315 

62.9208 

67 . 5656 

21 

0.9275 

62.5666 

67.6339 

22 

0,9183 

61.7857 

67.2827 

23 

0.9168 

61.6139 

67.2038 

29 

0.9136 

61.6339 

67.3677 

25 

0.9100 

61.6691 

67.5251 

26 

0.9126 

61.5941 

67.6066 

27 

0.9166 

61.9746 

67.6161 

2b 

0.9263 

62.6850 

67.6588 

29 

0.9339 

63.0337 

67.6981 

30 

0.9630 

63.5895 

67.6322 

31 

0.9516 

66.1703 

67.6322 

32 

0.9590 

66.7050 

67.6715 

33 

0.9703 

65.3221 

67.3202 

39 

0.9799 

65.8936 

67.2673 

35 

0.9886 

66.5602 

67.2929 

36 

0.9907 

56.9381 

67.5699 

37 

0.9930 

67.0170 

67.6901 

38 

0.9986 

67.6902 

67.5687 

39 

1.0000 

67.5682 

67.56*2* 

60 

1.0023 

67.7668 

67.6091 

61 

1.0000 

b7 • 561 6 

67.5414 

62 

0.9988 

67.5618 

67.6203 

63 

0.9983 

67.5021 

67.6198 

66 

1.0012 

67,9393 

67.8612 

65 

1.0023 

67.6257 

67.6686 

66 

1.0029 

67.6366 

67.6387 

67 

1.0006 

67.7502 

67.7110 

66 

1.0012 

67.7503 

67.6720 

69 

0.9983 

67.6813 

67.5992 

50 

1.0061 

67.8608 

67.5663 

51 

1.0017 

67,5395 

67.4215 

DRAG 

COtPP  IS 

0.0253 

BARON 

PRESS  IS 

28.960  IN.  HG 

SURVEY  N 

49 
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TABLE  D6f 

Survey  51,  a  =  10.36°. 


PSN 

RATIO 

VEL 

VRec 

l 

0.9953 

67.7351 

67.8489 

2 

1.0000 

6  7 . 74  7d 

67.7478 

3 

1.0000 

67.7160 

67.7160 

4 

1.000b 

67.5315 

67.4934 

5 

1.0017 

67.5060 

67.3914 

6 

1.0011 

67.6334 

67.557 2 

7 

C  .9994 

67.5635 

67.6016 

8 

1.0023 

67.6342 

67.5317 

9 

0.999*, 

67.4997 

67.5379 

10 

1  .0006 

67.6653 

67.6272  * 

11 

0.9955 

67.2890 

67.5948 

l  z 

0.9915 

67.1035 

67.6776 

13 

0.96ol 

66.8790 

67.6841 

14 

0.984b 

62.1710 

63.1439 

15 

0.97P8 

66.C304 

67.4606 

lb 

0.9704 

o5.8869 

67.8933 

17 

0.9594 

65.1049 

67.8617 

18 

0.9514 

64.3127 

67.6009 

19 

0.9427 

o  3 .° 1 6  2 

67.7982 

20 

0.9322 

63.0675 

67.6514 

21 

0.9195 

o2.2449 

67.6960 

22 

0.9131 

61.7069 

67.5812 

23 

0.9111 

61.5449 

67.5492 

24 

0.9085 

61.3346 

67.5110 

25 

0.9071 

61.3886 

67.6769 

26 

0.9058 

61.3043 

67.6769 

27 

0.9116 

61.7721 

67.7595 

28 

0.9101 

62.3138 

67.7977 

29 

0.9303 

62 .9446 

67.6575 

30 

0.9365 

63.4416 

67.7402 

31 

0.9460 

04.0161 

67.6704 

32 

0.9543 

64.5854 

67.6766 

33 

0.9643 

65.2292 

67.6445 

34 

0.9738 

66.0122 

67.7914 

35 

0.9789 

66.2592 

67.6892 

36 

0.9857 

66.6547 

67.6138 

37 

0.9926 

67.0544 

67.5547 

38 

0.9977 

67.4780 

67.6314 

39 

0.9977 

67.4459 

67.5994 

40 

1 .0000 

67.5165 

67.5165  * 

41 

1.0011 

67.7465 

67.6700 

42 

1.0000 

27.8397 

27.8397 

43 

1.0017 

67.8048 

67.6901 

44 

1.0011 

67.6947 

67.6130 

45 

1.0011 

67.8418 

67.7654 

46 

1.0017 

67.8479 

67.7332 

47 

1.0011 

67.7841 

67.7076 

48 

1.0000 

67.7138 

67.7138 

49 

1.0017 

67.7458 

67.6308 

50 

1.0011 

67.6433 

67.5665 

51 

1.0034 

67.7643 

67.5343 

Oft  AG 

COEFF  IS 

0.0280 

BAftOft 

PRESS  IS 

29.140  IN.  HG 

SURVEY  NO 

APPENDIX  E 
ERROR  ANALYSIS 
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Two  types  of  data  were  measured  during  the  course  of  the  investiga¬ 
tion:  (1)  data  measured  by  either  of  the  two  force  balances  and 

(2)  data  measured  by  the  wake  traverse  procedure.  In  the  force  balance 
phase  of  the  study,  there  were  several  intermediate  measurements  and 
calculations  involved.  The  air  density  was  determined  from 


p  =  0.002378(bar°-me^^f-adl-^)( 


519 


temperature  reading 


) 


(El) 


2 

Barometric  pressure  was  measured  to  ±0.01  in.  Hg  (±0.71  lb/ft  ). 
Temperature  was  measured  electronically  and  accurate  to  ±0.01°  F. 
Velocity  was  determined  from 


V  = 


(E2) 


where  AP  =  pressure  difference  across  the  transducer  and  p  «  density. 

2 

AP  was  measured  as  a  voltage  and  converted  to  lb/ft  .  Individual  volt¬ 
age  measurements  were  obtained  by  integrating  voltage  outputs  of  the 
voltmeter  over  a  10-sec  interval.  Voltage  measurements  of  AP  were  accu¬ 
rate  to  ±0.003  volt.  Forces  as  sensed  by  the  strain  gages  were  derived 
from 


FA,B<lb> 


(~rr)Volts 


A, B  volt 


A,B 


(E3) 


where  K.  _  *  calibration  constants  for  ribbons  A  and  B.  Calibration 
A,B 

constants  were  verified  by  using  the  balance  system  to  measure  known 
weights.  The  balance  system  was  able  to  measure  known  standard 
calibration  weights  to  within  ±0.6%  of  their  value.  The  strain  gages 
were  calibrated  numerous  times  with  a  resulting  variation  of 
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±0.000135  lb/volt.  The  airfoil  dimensions,  span  and  chord,  were  accu¬ 
rate  to  ±0.001  in.  Equations  (El),  (E2),  and  (E3)  were  combined  to  gen¬ 
erate  lift  and  drag  force  coefficients  in  the  form 


C 


force 


_ Force _ 

(l/2)pV^ (span) (chord) 


(E4) 


The  airfoil  angle  of  attack  was  measured  with  an  artillery  gunner's 
quadrant  in  mils  (1°  =  17.78  mils).  Angles  were  measured  with  an  accu¬ 
racy  of  ±0.01  mil  or  ±0.0056°. 

The  individual  accuracies  of  the  elements  making  up  Equation  (E4) 

can  be  used  to  examine  a  worst-case  condition  for  a  typical  C  or  C 

L  U 

calculation.  A  relative  error  for  each  of  the  terms  in  Equation  (E4) 
can  be  obtained  by  using  logarithmic  differentiation.  For  example. 
Equation  (El)  can  be  expressed  as 


P 


(El) 


where  B  <=  barometer  reading,  T  *  temperature  reading,  and  K  *  constant. 
Equation  (El)  can  also  be  expressed  as 


inp  ■  £nK  +  fcnB  -  fcnT  .  (El) 

Differentiating  both  sides  of  this  form  of  Equation  (El)  relates  the 
relative  errors  among  the  various  terms  in  Equation  (El) , 


Ap_  _  AB  AT 
p  "  B  “  T 


(Ela) 


The  sign  of  the  AT/T  term  can  be  +  or  Therefore,  for  a  maximum  error 
term,  Equation  (Ela)  becomes 
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Ap  AB  AT 

p  =  B  T 


(Elb) 


Performing  the  same  operation  on  Equations  (E2)  and  (E3) ,  one  gets 


AV  =  Ap  A  (AP) 
V  2p  2  AP 


(E2a) 


and 


AF,  n  AK,  _  AVolt  _ 

A,B  _  A,B  _ A,B 


A,B 


A,B 


Volt 


A,B 


(E3a) 


Performing  a  similar  operation  on  Equation  (E4)  and  substituting  (Elb) , 
(E2a) ,  and  (E3a) ,  one  gets  an  expression  for  the  overall  relative  error 
for  CL  or  CD, 


AK 


A,  B 


A,B 


AVOltA,B  2  \  +  2AT  Ab  Ac  A(_AP) 

Volt  _  ZCB  }  +  T  b  c  AP 

A  ,  D 


(EAa) 


where  b  =  span  and  c  =  chord. 

Inserting  the  specified  or  experimentally  determined  accuracy  data 
for  the  individual  components  of  Equation  (E4a)  on  a  term-by-term  basis, 
one  finds  the  following  worst-case  relative  error: 

ACf  0.0003  lb /volt  0.0003  volt  .0.02  in.  Hg  +  °-3°F  ) 

Cf  *  0.176  lb/volt  0.9  volt  ^30.0  in.  Hg;  530. 0°F 

0.002  in.  0.002  in.  0.01  pai 
18.375  in.  18.00  in.  1.00  psi 


0.0158 
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The  corresponding  absolute  error,  AC 


f* 


can  now  be  determined  from 


ACf  =  (^-i)Cf  .  (E5) 

For  a  typical  =  0.0250  corresponding  to  a  voltage  of  0.900  volts,  the 
cumulative  system  error  would  result  in  a  final  coefficient  error  of 

ACd  =  (0.0158)(0.025)  =  0.0004 


which  is  four  counts  of  drag. 

As  an  independent  assessment,  the  Student-t  distribution  was  con¬ 
sidered  to  appraise  statistically  the  possible  size  of  Cforce  errors. 
The  t-distribution  is  expressed  as 


t 


(E6) 


where  X  and  s  are  the  mean  and  standard  deviation  of  a  random  sample  of 
xize  n  from  a  population  which  has  a  mean  y,  a  standard  deviation  o,  and 
which  can  be  approximated  closely  by  a  normal  curve.  For  each  angle  of 
attack,  three  measurements  of  C  and  Cn  were  made.  Rearranging 

la  D 

Equation  (£6) ,  one  obtains 


ts 

& 


X  - 


V 


(E7) 


If  a  value  for  t  is  selected  from  a  t-distribution  table  corresponding 
to  a  probability  of  90%,  Equation  (E7)  will  give  a  corresponding  error. 
Then,  for  a  sample  of  n  measurements,  one  can  assert  with  a  90% 


135 


probability  that  the  mean  of  that  sample  has  an  error  of  less  than  that 
computed  by  Equation  (E7). 

For  example,  consider  data  measured  during  a  test  using  the  18-in. - 
chord  airfoil  without  disk  for  an  angle  of  attack  of  6.17°,  three  meas¬ 
urements  of  CD  were  made;  namely,  0.01582,  0.01553,  and  0.01577.  The 
sample  size  (n)  is  3.  The  sample  mean  (X)  is  0.01571.  The  standard 
deviation  (s)  is  given  by 


s  = 


z 


i  =  1 


(E8) 


For  this  example, 

» .  Jh^h  tlnMh  t  hm^h .  o.oooif,  . 

For  n  =  3  and  a  probability  of  90%,  a  t-distribution  table  (20)  yields 
t  =  2.920.  Substituting  these  values  in  Equation  (E7),  one  gets 

—  =  2- 920(0. 00016)  =  0.00027 

&  /3 

For  measurements  of  made  at  a  =  6.17°,  there  is  a  90%  chance  that  the 
mean  of  the  measurements  is  within  ±0.00027  of  the  correct  value.  As 
would  be  expected,  the  0.00027  error  is  within  the  bounds  of  the  worst- 
case  analysis.  This  type  of  analysis  was  performed  on  each  of  the 
various  force  balance  experiments.  For  a  given  experiment,  the  largest 
error  calculated  in  this  manner  was  selected  and  listed  at  the  head  of 
each  column  of  and  data.  See  the  tables  in  Appendix  D. 
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Drag  coefficients  obtained  from  the  wake  integration  process  are 
subject  to  the  previously  discussed  errors  associated  with  velocity 
measurement.  Since  the  airfoil  wake  is  measured  as  a  series  of  velocity 
ratios,  errors  in  velocity  tend  to  cancel  out.  There  is  an  error  asso¬ 
ciated  with  use  of  the  trapezoidal  rule  of  integration.  For  a  function 
of  y  =  F(x),  integrated  from  a  to  b  and  divided  into  trapezoidal  sub¬ 
intervals  of  width  h,  an  expression  for  the  so-called  global  error  is 
given  by  Reference  (2l)  as 

|E|  <.  -  (~-5-)h2f"(0  •  (E8) 


If  £  is  chosen  in  the  interval  (a,b)  such  that  it  is  a  maximum  value, 

then  Equation  (E8)  gives  the  upper  bound  of  the  error. 

Pope  (l3)  suggests  that  the  curve  of  an  airfoil  wake  profile  may  be 

2 

approximated  by  a  (cosine)  function.  Consider  Survey  No.  45 
(Appendix  D) ,  where  (b  -  a)  =  2.0  in.  (0.1667  ft)  and  h  «  0.1  in. 

(0.0083  ft).  The  wake  profile  can  be  expressed  as 


cos2((^4y)x)  - 

n  ' v  n  n  ma V  V  L/ 2  ' 


qo°X 


q  max 
o 


(E9) 


where  L  »  (b  -  a)  *  width  of  the  wake  (integration  limits). 

The  quantity  in  the  first  set  of  brackets  represents  the  maximum 
amplitude  of  the  wake  profile  curve.  Equation  (E9)  models  a  curve 
similar  to  that  shown  on  Figure  14.  Substituting  (E9)  into  (E8)  gives 
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To:  Distribution  List 

Subject:  Forwarding  of  Errata  for  ARL  TM  80-44 


Please  replace  Page  82  of  ARL  TM  80-44,  "A  Method  of  Correcting 
for  the  Effects  of  the  Sidewall  Boundary  Layer  in  Two-Dimensional 
Airfoil  Testing,"  by  Paul  P.  Jacobs,  Jr.,  dated  March  31,  1980. 

An  error  in  definition  was  found  in  Paragraph  1  and  also  a 
typographical  error  was  made  in  Equation  (13). 
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where  C£  ■  dC^/da  and  C£  ■  (dC£/ dct)Q.  As  the  experimental  data  show, 

*  “o 

the  drag  correction  must  be  zero  at  C#  -  0.  Thus,  the  C,  contribution 

should  be  directly  proportional  to  the  AC.  term.  The  assumed  C.  cor- 

a  a 

rection  can  be  the  term  which  "Individualizes"  the  correction  to  spe¬ 
cific  airfoils  if  it  is  set  equal  to  Cd  (the  drag  coefficient  at  zero 

o 

lift)  of  the  tested  airfoil.  The  correction  now  takes  the  form 

cdVc»/ci  >1/2 

o  a  a 

ACd  -  K - 2 -  f  (Alt)  .  (12) 

((fi/c)/0.l)(t/c)4 

At  this  point.  Equation  (12)  can  be  applied  to  the  data  of 
Figures  40-42.  On  a  trlal-and-error  basis,  setting  K  ■  0.00015  and 
f(AR)  -  1/ /ar  ,  one  finds  for  the  final  correction  that 

0.00015(cd  C£(C£  /C£  )1/2) 
o  a  a 

ACd - 2 -  .  (13) 

((6/c)/0.l)(t/c)4vfiR 

Shown  in  Figures  43-45  are  the  results  of  the  three  experiments 
with  Equation  (13)  applied  to  the  data.  The  correction  works  extremely 
well  for  the  ARL/PSU  data.  There  is  also  exceptionally  good  agreement 
when  it  is  applied  to  Kermeen's  data.  With  Kermeen's  data,  there  is  a 
slight  deviation  at  the  higher  C£  values.  It  must  be  noted  that  In  the 
range  C£  •  1.1  Kermeen's  curve  exceeds  the  published  NACA  data  by 
approximately  90Z.  Another  consideration  is  the  advances  made  in 
instrumentation  between  Kermeen's  tests  and  the  present  study.  Also 
noteworthy  is  the  fact  that  the  Reynolds  numbers  of  the  two  curves 


